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MASSACHUSETT 


ABSTRACT 


This  report  summarizes  the  current  level  of  understanding  of  the  performance  of  arrays 
of  real  radiators.  The  extent  and  nature  of  theeffects  ofmutual  coupling  (element  pat¬ 
tern  distortion,  element  impedance  variation  with  scan  angle,  and  polarization  varia¬ 
tion  with  scan  angle)  in  array  antennas  are  investigated  for  a  variety  of  array  geometries 
and  for  several  types  of  radiating  elements.  Both  finite  and  infinite  arrays  of  regu¬ 
larly  spaced,  uniformly  illuminated,  and  progressively  phased  elements  are  consid¬ 
ered.  The  effects  of  coupling  on  unequally  spaced  arrays,  arrays  with  coupled  feed 
networks,  and  circularly  polarized  arrays  are  also  discussed.  Finally,  the  results  of 
a  study  of  the  effects  of  coupling  on  the  radiation  patterns  of  multiple-beam  optical- 
type  antennas  are  presented.  Although  most  of  the  numerical  results  are  based  on 
thin,  dipole  radiating  elements  with  and  without  ground  planes  (the  dipole  without  a 
ground  plane  is  the  exact  dual  of  an  arrayof  slots  in  a  ground  plane),  both  theoretical 
and  experimental  investigations  of  the  sensitivity  of  the  results  to  the  element  type  are 
included.  Emphasis  is  placed  on  the  results  of  various  analyses  and  their  implications 
to  the  array  designer.  The  derivations  of  some  of  the  more  important  results  are  briefly 
outlined,  but  in  most  instances,  only  the  relevant  conclusions  are  given. 
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FOREWORD 


The  primary  objective  of  this  report  is  to  provide  the  working  engineer  with  a  knowledge 
of  the  effects  of  mutual  coupling  on  array  performance  —  what  the  effects  are  and  how  serious 
they  are.  The  theory  and  the  necessary  analytical  tools  for  reproducing  and  extending  the 
results  are  outlined  and  referenced,  but  the  emphasis  is  on  results;  i.e.,  to  what  extent  do 
the  properties  of  an  array  depend  on  the  kind  of  radiating  elements  used  and  the  way  they 
are  excited? 

The  report  thus  has  shortcomings  that  would  be  of  interest  to  the  well-rounded  engineer:  it 

ignores  certain  very  fundamental  results*  that  have  not  yet  proven  of  practical  value,  and  it 
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ignores  entirely  the  question  of  how  to  measure  experimentally  the  effects  discussed.  Al¬ 
though  the  purpose  of  this  report  is  tutorial,  many  of  the  results  are  appearing  here  for  the 
first  time. 

The  origins  of  this  report  can  be  traced  to  1961  when  one  of  the  authors  prepared  notes  on 
array  antenna  theory  for  a  summer  course  in  "Elements  of  High  Powered  Radar  Design"  at 
M.I.T.  In  1963,  the  revised  notes  were  issued  as  a  treatise  on  array  theory. ^  At  that  time, 
the  author  was  dissatisfied  with  the  state  of  knowledge  on  mutual  coupling;  subsequently,  both 
authors  of  the  present  report  have  devoted  a  great  deal  of  time  trying  to  improve  the  under¬ 
standing  of  such  phenomena.  They  were  also  fortunate  to  establish  and  maintain  person-to- 
person  communication  with  many  others  who  have  become  interested  in  this  area  and  whose 
names  appear  frequently  in  the  references.  As  a  result  of  these  efforts,  significant  progress 
was  made  which  thoroughly  outdated  the  chapter  onmutual  coupling  in  Ref.  6  in  just  a  few  years. 

The  motivation  for  this  updating  of  the  authors'  writings  on  mutual  coupling  research  is  due  to 
Professor  Curt  Levis,  Director  of  the  Ohio  State  University  Antenna  Laboratory.  At  his  invi¬ 
tation,  the  authors  gave  lectures  during  the  summers  of  1965  and  1966  as  part  of  the  Univer¬ 
sity's  curriculum  on  "Recent  Advances  in  Antenna  and  Scattering  Theory."  The  1965  course 
notes  formed  the  original  draft  for  this  report.  For  the  1966  course,  the  notes  underwent  a 
major  revision.  Since  that  time,  only  a  few  corrections  have  been  made. 

The  authors  are  indeb.ed  to  Professor  Levis  for  the  inspiration  to  write  this  report  and  for 
permission  to  publish  it  as  a  Lincoln  Laboratory  Technical  Report.  They  are  also  apprecia¬ 
tive  of  those  who  supported  the  work  at  Lincoln  Laboratory,  namely,  the  U.  S.  Air  Force 
and,  more  recently,  the  Advanced  Research  Projects  Agency.  The  authors  especially  thank 
Lt.  Colonel  John  C.  Toomay  who  was  program  director  atARPAandwho  has  encouraged  broad 
research  of  this  type,  insisting  that  the  results  be  written  in  a  manner  intelligible  to  those  not 
specializing  in  the  field. 

It  would  be  impossible  to  list  all  those  whose  research  contributed  to  the  body  of  knowledge 
described  in  this  report.  The  references  do  this  to  some  extent,  but  do  not  give  credit  to  the 
engineer  whose  experimental  work  is  not  published  in  the  open  literature.  In  this  category,  the 
authors  especially  acknowledge  the  contributions  of  L.  Schwartzmann  of  Sperry  Gyroscope 
Company  and  R.  Tang  of  Hughes  Aircraft  Company  whose  observations  and  experiments  have 
been  helpful  in  verifying  and  extending  the  theory  of  mutual  coupling  effects. 

For  the  interested  reader,  the  authors  recommend  Microwave  Scanning  Antennas  —  Vol,  II:  Ar¬ 
ray  Theory  and  Practice,  R.  C.  Hansen,  editor  (Academic  Press,  New  York,  1966).  Chapters 
2  through  4  by  A.  A.  Oliner  and  R.G.  Malech  cover  much  the  same  material  as  this  report  but 
from  a  somewhat  different  viewpoint  in  some  cases. 
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MUTUAL  COUPLING  IN  ARRAY  ANTENNAS 


I.  INTRODUCTION 

A.  Qualitative  Explanation  of  Coupling  Effects 

It  is  intrinsic  to  the  nature  of  antennas  that  when  two  antennas  are  in  proximity  and  one  is 
transmitting,  the  second  will  receive  some  of  the  transmitted  energy,  with  the  amount  depend¬ 
ent  on  their  separation  and  relative  orientation.  Even  if  both  antennas  are  transmitting,  they 
will  simultaneously  receive  part  of  each  other's  transmitted  energy.  Furthermore,  antennas 
rescatter  a  portion  of  any  incident  wave  and  thus  act  like  small  transmitters  even  when  they  are 
nominally  only  receiving.  The  result  is  that  energy  interchange  between  a  particular  element 
of  an  array  and  a  remote  point  occurs  not  only  by  the  direct  path,  but  also  indirectly  via  scat¬ 
tering  from  the  other  antennas  of  the  array.  This  effect  is  a  manifestation  of  the  "mutual  cou¬ 
pling"  that  exists  between  array  antennas.  It  is  not  usually  a  negligible  effect  and  complicates 
the  design  of  such  antennas. 

Before  attempting  a  quantitative  examination  of  coupling  effects,  let  us  elaborate  on  the 
foregoing  physical  picture  by  examining  (1)  elements  of  a  transmitting  array  and  (2)  elements 
of  a  receiving  array  (which  may,  in  fact,  be  the  same  array  at  a  different  time). 

i.  Coupling  in  Transmitting  Arrays 

Imagine  the  two  antennas  shown  in  Fig.  1  as  being  two  of  many  in  an  array.  The  generator 
attached  to  antenna  n  sets  up  a  wave  traveling  from  the  generator  toward  the  antenna  indicated 
by  the  arrow  labeled  (0)  in  the  figure.  Part  of  this  energy  is  radiated  directly  into  space  (i), 
while  a  part  is  coupled  to  the  other  antennas  of  the  array  —  in  this  particular  case,  to  antenna  m. 
The  field  incident  upon  antenna  m  causes  current  flow  in  that  antenna  which  reradiates  some  of 
the  received  energy  (3)  and  also  launches  a  wave  toward  the  generator  of  antenna  m  (4).  Of  the 
energy  rescattered  (3),  some  is  reradiated  directly  into  space  and  some  in  turn  couples  again  to 
other  elements,  and  so  forth.  If  antenna  m  is  also  being  excited  by  its  own  generator,  the  en¬ 
ergy  rescattered  from  antenna  m  due  to  the  generator  excitation  of  antenna  n  adds  vectorially 
to  the  energy  from  generator  m,  altering  the  amplitude  and  phase  of  apparent  excitation  of  an¬ 
tenna  m  in  a  manner  dependent  on  the  output  of  generator  n. 

Thus,  the  total  contribution  to  the  far-field  pattern  of  a  particular  element 
in  the  array  depends  not  only  on  the  excitation  furnished  by  its  own  generator 
(the  direct  excitation),  but  also  on  the  total  parasitic  excitation,  which  de¬ 
pends  on  the  couplings  from  and  the  excitation  of  the  other  generators. 

The  wave  which  is  traveling  toward  the  generator  (4)  of  antenna  m  adds  to  any  reflection 
from  that  antenna  because  of  mismatch  between  the  generator  and  the  antenna.  For  the  case  of 
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Fig.  ).  Coupling  paths  from  antenna  n 
to  antenna  m  when  transmitting. 


Fig.  2.  Coupling  paths  from  antenna  m 
to  antenna  n  when  receiving. 
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principal  interest  here,  in  which  the  element  excitations  are  coherent,  the  wave  (4)  due  to  gen¬ 
erator  n  differs  from  the  reflected  wave  in  antenna  m  because  of  its  own  generator  mismatch 
only  in  phase  and  amplitude;  exactly  how  depends  on  the  coupling  between  the  two  antennas  and 
the  complex  amplitude  of  the  output  of  generator  n.  The  net  effect  of  the  vector  sum  of  these 
two  waves  to  an  observer  looking  at  antenna  m  from  its  generator  is  the  same  as  if  the  imped¬ 
ance  of  antenna  m  were  changing  as  a  function  of  the  excitation  and  location  of  antenna  .41.  In 
fact,  it  is  common  pracvice  to  model  this  effect  of  coupling  on  such  arrays  as  a  change  in  the  ap¬ 
parent  driving  impedance  of  the  elements  (this  change  is  often  termed  the  "mutual  impedance 
variation" ). 

While  not  entirely  without  its  confusing  aspects,  the  concept  of  a  changing  driving  imped¬ 
ance  is  useful.  For  example,  for  a  particular  set  of  element  excitations  and  locations  in  an  ar¬ 
ray,  the  generator  impedance  that  is  optimum  for  each  element  (in  the  sense  of  maximizing  the 
radiated  power)  would  be  a  "match"  to  an  impedance  at  the  antenna  terminals  that  would  set  up  a 
reflected  wave  numerically  equal  in  phase  and  amplitude  to  the  backward  traveling  waves  induced 
because  of  the  coupling.  Such  a  generator  achieves  maximum  power  transfer  in  an  array,  even 
though  it  is  not  the  impedance  which  is  a  match  to  what  one  would  measure  looking  into  a  single 
antenna  without  the  other  antennas  excited. ^  Since  the  coupled  waves  depend  on  the  excitation  as 
well  as  the  placement  of  the  other  antennas,  the  impedance  selected  is  optimum  only  for  that  set 
of  conditions.  Thus  we  arrive  at  one  of  the  principal  aggravations  of  coupling  in  arrays  whose 
element  excitations  are  varied  (such  as  in  electronic  scanning  arrays): 

The  generator  impedance  that  would  maximize  the  array  radiation  efficiency 
(gain)  varies  with  the  array  excitation. 

2.  Coupling  in  Kcceiving  Arrays 

Figure  2  illustrates  the  coupling  paths  for  a  receiving  array  in  which  the  receivers  are  rep¬ 
resented  as  passive  loads.  Again,  the  figure  depicts  two  elements  of  a  large  number  assumed 
to  be  present  in  an  array. 

Assume  that  there  is  a  plane  wave  incident  as  indicated  by  the  arrows  labeled  (0)  from  a 
direction  such  that  the  incident  wave  will  strike  antenna  m  first.  The  field  incident  on  antenna 
m  causes  current  to  flow  on  this  antenna  which  launches  a  wave  into  its  feed  (1)  and  also  rescat¬ 
ters  some  of  the  energy  into  space  (2)  and  into  adjacent  antennas  (3).  The  rescattered  wave  from 
antenna  m  (3)  adds  vectorially  at  antenna  n  to  the  wave  directly  incident  from  space. 


t  In  order  to  minimize  confusion,  we  adopt  the  following  terminology: 

(a)  Antenna  impedance:  the  impedance  looking  into  a  single  isolated  element. 

(b)  Passive  driving  impedance:  the  impedance  looking  into  an  element  of  an 
array  with  all  other  elements  in  place  and  passively  terminated  (in  their 
normal  generator  internal  impedances,  unless  otherwise  specified). 

(c)  Active  driving  impedance:  the  impedance  looking  into  an  element  of  an 
array  with  all  other  antennas  in  ploce  and  excited  (excitation  must  be 
specified). 

Since  impedance  (b)  is  of  minor  importance  and  in  most  practical  cases  differs  only  slightly  from  (a),  the  term 
"driving  impedance"  alone  will  be  understood  to  refer  to  the  octive  driving  impedance. 


Thus  the  total  input  energy  to  each  antenna  from  the  incident  wave  is  the 
vector  sum  of  the  waves  coming  in  directly  from  space  as  well  as  those 
coupled  parasitically  from  the  other  antennas:  it  is  dependent  on  the  rel¬ 
ative  location  of  all  the  elements  in  the  array. 

Since  the  relative  amount  of  energy  absorbed  and  energy  reradiated  from  any  antenna  de¬ 
pends  on  its  match  to  its  terminating  impedance,  the  total  input  energy  to  each  antenna  depends 
on  the  terminating  impedances  of  all  the  other  antennas.  In  fact,  for  maximum  extraction  of 
energy  from  the  passing  wave,  we  would  like  to  choose  the  terminating  impedance  of  the  elements 
in  such  a  manner  that  we  minimize  the  total  energy  backscattered  into  space  (2).  Consequently, 
we  should  actually  mismatch  the  receiver  relative  to  the  ar.tenna  impedance  in  order  to  set  up  a 
reflection  from  the  receiver  back  to  the  antenna  (4)  to  cancel  the  wave  that  would  have  been  re¬ 
scattered  had  the  receiver  been  matched  to  the  actual  impedance  of  each  antenna: 

The  proper  impedance  depends  on  the  placement  and  excitation  of  the  other 
elements.  In  fact,  the  optimum  receiver  input  impedance  is  precisely  the 
same  as  the  optimum  generator  impedance  for  the  same  array  for  trans¬ 
mitting  energy  in  the  same  direction. 

This  fact  follows  directly  from  the  principle  of  reciprocity.  It  can  be  inferred  directly,  although 
it  is  by  no  means  a  simple  exercise. 

It  should  now  be  evident  that  mutual  coupling  plays  an  important  role  in  the  performance  of 
array  antennas.  It  may  not  be  so  evident  that  understanding  and  coping  with  its  effects  can  be 
reduced  to  a  reasonable  problem.  The  remainder  of  this  section  is  devoted  to  a  detailed  anal¬ 
ysis  of  the  ways  in  which  mutual  coupling  depends  on  the  design  parameters  of  the  array  and  its 
elements. 

B.  Effects  of  Mutual  Coupling  on  Array  Performance 

The  exact  extent  and  nature  of  the  effects  of  coupling  on  array  performance  depend  on  (1) 
the  type  of  antenna  and  its  design  parameters,  (2)  the  relative  placement  of  the  elements  in  the 
array,  (3)  the  type  of  feed  used  to  excite  the  elements  and  the  design  parameters  thereof,  and 
(4)  the  range  of  relative  excitations  employed  (the  scan  volume  of  the  array).  In  the  following 
subsections,  we  will  examine  quantitatively  the  extent  of  these  effects  and  their  dependence  on 
the  factors  mentioned.  Before  doing  so,  however,  let  us  point  out  the  practical  consequences 
of  tnese  effects.  These  fall  into  two  classes:  (t)  those  effects  which  arise  from  the  apparent 
variation  in  element  driving  impedance  and  (2)  those  which  arise  from  the  multipath  nature  of 
the  route  followed  by  the  energy  from  each  generator  to  the  far  field  or  from  the  incident  wave 
to  each  receiver.  These  cause  distortion  of  the  array  far-field  pattern. 

1.  Impedance  Variation  and  Its  Subsidiary  Effects  on  Equipment 

The  apparent  variation  in  element  impedance  generally  leads  to  a  variation  in  array  effi¬ 
ciency,  since  it  is  a  practical  impossibility  to  match  the  element  driving  impedance  for  all  con¬ 
ditions  of  excitation.  In  addition,  the  apparent  mismatch  produces  precisely  the  same  effects  as 
a  real  mismatch  on  such  auxiliary  equipment  as  receivers,  transmitters,  and  transmission  lines. 
Specifically,  it  can  contribute  to  increasing  the  noise  figure  of  the  receivers  and  distort  the  re¬ 
ceiver  transfer  coefficient  in  both  phase  and  amplitude  (leading  in  severe  cases  to  receiver  in¬ 
stability).  On  transmission,  the  transmitter  transfer  characteristic  :an  be  similarly  altered. 
Also,  on  transmission,  there  are  the  possibilities  of  voltage  breakdown  and  overheating  in  the 
transmitter  output  lines  due  to  standing  waves,  and  degradec  transmitter  output  due  to  detuning. 
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2.  Pattern  Variation  Effects 

The  multipath  route  the  energy  follows  because  of  mutual  coupling  causes  the  patterns  of 
arrays  to  differ  from  those  that  would  be  predicted  on  the  basis  of  simple  theory  using  noninter¬ 
acting  elements.  In  very  large  regular  arrays  (those  in  which  the  elements  are  placed  at  regular 
intervals  on  a  grid  and  which  have  a  sufficient  number  of  elements  so  that  edge  effects  can  be  ig¬ 
nored),  the  usual  pattern  distortion  effect  is  a  simple  scaling  up  and  down  of  the  relative  amplitude 
of  the  patterns  while  preserving  the  relative  pattern  shape.  However,  for  irregular  arrays  or  for 
small  regular  arrays  with  dominating  edge  effects,  the  relative  pattern  differs  substantially,  in 
ways  which  are  often  quite  complicated,  from  the  pattern  computed  ignoring  coupling  effects. 

In  antennas  with  elements  that  can  support  more  than  one  orthogonal  sense  of  polarization, 
it  is  possible  for  the  coupling  to  excite  the  polarization  sense  which  is  not  directly  excited  by  the 
generators,  thus  causing  depolarization  of  the  signal. 

Both  the  change  in  polarization  and  the  change  in  efficiency  due  to  mismatch  effects  cause 
the  gain  of  the  array  to  change  as  a  function  of  the  relative  excitation  of  the  elements;  for  ex¬ 
ample,  the  gain  varies  with  scan  angle  (element -to-element  phasing)  in  phase-scanned  arrays. 

n.  COUPLING  IN  INFINITE  REGULAR  ARRAYS 

As  an  aid  to  analysis  and  understanding,  it  is  helpful  to  analyze  an  idealized  array  which, 
although  a  strict  physical  impossibility,  is  close  enough  to  the  real  world  to  provide  directly  use¬ 
ful  answers  in  many  cases  of  practical  importance,  and  insight  into  others. 

A.  Infinite  Array  Concept  and  Its  Utility 

The  most  useful  idealized  array  model  for  coupling  studies  is  usually  referred  to  simply  as 
an  "infinite  array."  This  means  an  array  with  (1)  all  elements  placed  at  regular  intervals  either 
in  a  straight  line  or  on  a  flat  surface  of  infinite  extent,  (2)  all  elements  identical,  (3)  the  ampli¬ 
tude  of  element  excitation  equal  from  one  element  to  another  (uniform  amplitude  illumination), 
and  (4)  the  relative  phasings  of  the  elements  differing  at  most  by  a  term  which  is  linear  in  two 
orthogonal  directions  across  the  array.  Mathematically,  if  the  array  lies  in  the  x-y  plane  of  a 
coordinate  system  such  as  that  of  Fig.  3  (centering  one  element  for  convenience),  the  elements 


Fig.  3.  Generalized  planar  array  geometry. 
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t 

are  located  at  the  intersections  of  a  rectangular  grid  with  spacings  Dx  and  D^.  Thus,  it  is  con¬ 
venient  to  represent  elements  by  a  double  set  of  indices,  so  that  the  mn**1  element  is  located  at 

x  =  mD  ,  y  =  nD  .  The  element  excitations  differ  only  by  a  phase  term 
x  y 

*mn  =  m“  + 

where  a  and  0  are  the  row  and  column  phase  increments. 

The  utility  of  this  model  is  in  the  simplifications  that  it  permits  in  mathematically  modeling 
and  analyzing  the  coupling  in  an  array.  As  will  be  seen,  the  behavior  of  such  an  infinite  array 
describes  fairly  accurately  the  behavior  of  most  of  the  elements  of  modest  to  large  arrays  that 
are  either  on  flat  surfaces  or  shallowly  curved  surfaces  with  a  smoothly  varying  amplitude  and 
phase  taper. 

The  nature  of  the  behavior  of  the  element  driving  impedance  as  a  function  of  scan  angle  is 
easily  assessed  in  an  infinite  array.  For  an  antenna  which  is  fed  by  a  single-mode  transmission 
line,  we  can  write  the  terminal  voltage  of  any  one  antenna  in  terms  of  the  current  flowing  in  the 
others  by  simultaneous  equations  of  the  form 


V  =  )  Z  I  (1) 

mn  mn,  pq  pq  '  ' 

where  the  equation  defines  Z_  as  the  ratio  of  the  terminal  voltage  at  antenna  mn  due  to  a 

mn,  pq 

unity  current  flowing  in  antenna  pq  when  all  other  antenna  currents  are  zero.  Hence,  the  Z„ 

mn,  pq 

are  termed  "mutual  impedances"  when  the  indices  mn  and  pq  are  different.  The  driving  imped¬ 
ance  of  the  mnth  antenna  is,  by  definition. 


nm 


nm 


nm 


=  Z  Z  znm,  pq 


m  n 


nm 


(2) 

(3) 


But  the  regularity  of  the  array  and  its  excitation  requires  that 


pq 


*00* 


j(pa+q/3) 


If  we  consider  the  central  element  for  notational  simplicity,  we  have 

ZD(“.fl)=  ZZZoo,p,ei(P“*q'J1 


m  n 


(4) 

(5) 


which  makes  apparent  the  fact  that  the  driving  impedance  for  such  an  array  is  simply  the  vector 
sum  of  ZQo  QO  and  the  so-called  "phased  mutual  impedances"  resulting  from  the  excitation  of  the 
other  antennas.  The  term  ZOQ  QQ  is  the  impedance  looking  into  the  central  element  with  all  others 
open -circuited  so  that  the  current  at  their  feed  points  is  zero.  If  so  terminating  the  antennas 
causes  the  current  to  be  zero  everywhere  on  the  element  (as  would  be  true  for  thin,  gap-fed  di¬ 
poles  and  slots,  if  the  reference  point  is  effectively  at  the  gap),  then  setting  I  =  0  is  physically 
equivalent  to  removing  the  pq  element,  and  Zqq  qq  is  numerically  the  impedance  of  a  single  iso¬ 
lated  element.  In  most  practical  cases,  it  is  very  nearly  the  same.  We  will  investigate  the  be¬ 
havior  of  the  mutual  impedances  further  in  Sec.  III. 


fThe  practical  and  important  cat*  of  a  triangular  grid  can  b«  analyzed  by  vacating  avary  other  site  of  a  rec¬ 
tangular  grid. 


Since  Eq.  (5)  was  obtained  without  reference  to  the  feed  system  supplying  the  element  excita¬ 
tions,  we  can  state  categorically  that  in  any  infinite  planar  array  of  single-mode  elements,  the 
driving  impedance  is  given  by  a  two-dimensional  Fourier  series  in  the  phasing  constants.  The 
coefficients  of  the  series  are  the  complex  mutual  impedances  between  elements.  The  ramifica¬ 
tions  of  this  result  will  be  explored  more  fully  below. 

A  second  virtue  of  the  infinite  array  model  is  that  its  periodic  nature  allows  us  to  i...ack 
boundary  value  problems  by  use  of  Fourier  series  (see  for  example  the  discussion  of  the  grating- 
lobe  series  in  Sec.  II-D). 

Finally,  since  the  impedance  behavior  of  an  antenna  in  an  infinite  array  of  single-mode  ele¬ 
ments  does  not  depend  on  the  details  of  the  circuit  feeding  the  array,  we  can,  by  this  artifice,  sep¬ 
arate  the  discussion  of  the  intrinsic  properties  of  antenna  elements  from  the  more  complex  prob¬ 
lems  of  the  interaction  between  the  elements  and  the  feeds.  For  a  finite  or  irregular  array,  we 
cannot  in  general  discuss  antenna  impedance  as  an  isolated  factor. 

B.  Other  Useful  Idealizations 

Several  other  idealizations  which  are  commonly  used  are:  (i)  The  concept  of  a  "single-mode 
antenna  element":  an  element  whose  boundary  conditions  are  such  that  any  current  that  flows  has 
to  conform  to  a  single  unique  spatial  distribution.  The  infinitesimally  thin  half-wavelength  dipole 
and  the  thin  resonant  slot  are  classical  examples  of  this  type  of  element.  (An  open  waveguide 
which  has  substantial  excitation  of  higher  order  modes  at  the  mouth  is  not  of  this  type.)  (2)  The 
concept  of  an  "invisible  feed."  This  idealized  feed  does  not  perturb  the  current  distribution  on 
the  antenna  and  has  no  coupling  to  the  other  antennas  nor  to  other  feeds,  nor  does  it  cause  arty 
rescattering.  (3)  If  the  elements  are  above  a  ground  plane,  the  ground  plane  is  assumed  to  be 
infinite  in  extent  and  to  have  infinite  conductivity. 


C.  Resistive  Sheet  Problem 


The  fundamental  nature  of  the  mutual  coupling  problem  (or,  more  precisely,  one  of  its  ram¬ 
ifications  -  the  change  ir.  impedance  of  an  array  with  scan  angle)  is  illustrated  by  examining  the 
behavior  of  the  reflection  coefficient  of  a  thin  resistive  sheet  of  infinite  extent  backed  by  an  open 
circuit7  so  that  there  are  no  fields  transmitted  through  the  sheet  [see  Fig.  4(a)].^ 

The  boundary  condition  imposed  by  such  a  sheet  is  that  for  z  =  0,  for  all  x  and  y,  the  ratio 
of  the  total  tangential  electric  field  to  the  total  tangential  magnetic  field  be  numerically  equal  to 
the  surface  impedance  of  the  sheet.  We  choose  the  surface  resistivity  of  the  sheet  to  be  that  of 
free  space:  {  =  iZOxohms,  so  that  the  sheet  is  an  exact  match  for  a  normally  incident  wave. 

If  we  now  move  the  source  through  an  angle  ©E  from  the  normal  in  the  y-z  plane  of  Fig,  4(b), 
we  can  express  the  spatial  dependence  of  the  incident  and  reflected  plane  waves  as 


E  1  -  exp  [-jk(y  sin  ©E  +  z  cos  ©E)J 


where  2  *  has  a  tangential  component 


E<!tang  =  Eo  co80E 


T  A*  pointed  out  in  Ref.  7,  the  open  circuit  U  not  realizable,  nevertheless  (he  model  is  useful  as  we  shall  see. 


7 


Fig.  4.  Resistive  sheet  reflection  geometry: 
(o)  incident  wove  situations,  (b)  E-plane  scan 
reflection,  (c)  H-pfane  scan  reflection. 
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cos©E  expl-jky  sin©E)  +  cos  ©f  expt-jky  sin©r) 


H<j  exp [—jky  sln©E)  -  expl-jky  sin©r] 


=  { 


For  each  wave,  we  also  have 

,i 


Hi_fo 

Ho-  i 


and  H  Z=-f- 


giving 


E^  (cos  ©E  -  1)  exp  [-jky  sin  ©E)  +  (cos  ©r  +  1)  exp  [—jky  sin©r]  =  0 


This  equation  can  obviously  have  a  solution  independent  of  y  only  if  ©E  =  ,©r-  Making  this  sub¬ 
stitution  and  solving  for  the  voltage  reflection  coefficient  gives 


r 


E  = 


1  -  cos  0_ 
_ EL 

1  +  COS  ©E 


(6) 


or 


(7) 


By  analogy  to  the  transmission-line  formula  for  the  normalized  load  impedance  r  in  terms  of 
the  reflection  coefficient 


r  = 


r  -  1 
r  +  1 


the  sheet  represents  an  impedance  to  the  wave  that  varies  as 

r  =  -* 

E  cos  ©E 

For  scan  in  the  H-plane,  the  assumed  direction  of  E-field  of  Fig.  4(c)  leads  to 


cos  -  1 

r  s  -  i J u,. 

H  cos  ©H  +  1 

(9) 

(10) 

rH  *  CO,SH  • 

(H) 

Thus,  even  for  such  a  simple  situation  as  a  wave  impinging  on  a  resistive  sheet,  there  is  a 
change  in  the  apparent  impedance  that  the  sheet  presents  to  the  wave  as  a  function  of  scan  angle. 

The  practical  importance  of  this  result  becomes  apparent  if  we  assume  that  the  resistive 
sheet  represents  the  limiting  case  of  an  array  of  infinitesimal  dipoles  with  infinitesimal  inter¬ 
element  spacings.  It  is  logical  to  expect  on  this  basis  that  the  reflection  coefficient  behavior 
represented  by  the  sheet  would  be  an  extreme  case.  It  might  also  seem  logical  that  extremely 
dense  packing  of  elements  would  represent  maximum  Impedance  variation;  however,  comparison 


9 


indicates  the  opposite  is  true.  The  resistive  sheet  impedance  variation  appears  to  represent  a 
lower  bound  on  impedance  variation  with  scan. 


0  >0  40  to  to 


ANGLE  (dt«) 


Fig.  5.  Reflection  coefficient  magnitude  vs  principal  plane  scan  angle. 

For  example.  Fig.  5  compares  the  computed  variation  in  reflection  coefficient  magnitude 
for  arrays  of  half-wavelength  dipoles  for  various  spacings  D  and  various  heights  s  above  a 
ground  plane  with  the  reflection  coefficient  behavior  of  the  resistive  sheet.^  (The  curves  for 
wider  spacings  are  discontinued  at  the  angles  where  major  lobes  of  the  array  are  equally  dis¬ 
posed  about  broadside;  the  behavior  just  retraces  the  curve  from  that  point.)  Note  that  the  re¬ 
flection  coefficient  for  a  given  scan  angle  in  the  principal  plane  becomes  monotonically  better 
as  the  elements  are  brought  closer  together.  As  a  confirmation,  computations  of  the  impedance 
behavior  for  a  central  element  of  a  large  array  of  short  dipoles  without  a  ground  plane  for  Dx  = 
Dy  =  0.1X  are  also  shown.  It  is  apparent  that  the  approximation  to  the  impedance  variation  of 
the  resistive  sheet  is  quite  good. 

The  second  important  consideration  is  the  base  from  which  this  change  takes  place;  in  the 
figure,  this  is  the  impedance  which  matches  the  array  at  broadside.  What  the  figure  does  not 
show  is  that  for  the  very  closely  spaced  short  dipoles,  the  absolute  impedance  -  particularly 
the  reactive  part  —  has  become  exceedingly  large,  although  the  change  in  impedance  is  quite 
small.  All  the  data  presented  assume  the  generator  impedance  is  matched  to  the  element  driving 
impedance  when  the  array  is  phased  for  broadside  reception  or  transmission.  To  do  so,  it  is 
necessary  to  tune  out  the  broadside  driving  reactance.  If  this  reactance  is  large  compared  to 
the  resistive  part,  the  tuning  is  extremely  critical;  such  is  the  case  for  short  dipoles  stacked 
close  together. 


t  Th*  height  s  in  Fig.  5  Is  th*  on*  that  minimizes  th*  mismatch  in  th*  principal  planes  of  scan. 
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D.  Grating -Lobe  Series  and  Its  Significance 

The  periodic  nature  of  the  infinite  array  permits  a  very  general  analysis  of  the  impedance- 
vs-scan  behavior  of  various  elements.  By  expressing  the  pov.er  radiated  and  the  power  stored 
in  the  vicinity  of  the  array  in  terms  of  the  current  or  field  distribution  on  a  typical  element,  the 
complete  impedance  behavior  as  a  function  of  inter-element  phasing  can  be  described.  The  re¬ 
sult,  although  in  the  form  of  an  infinite  series,  is  not  only  a  useful  computational  tool  for  ele¬ 
ments  with  known  current  distributions,  but  also  furnishes  insight  into  the  properties  of  elements 
that  affect  their  impedance  behavior. 

8  9 

Pioneering  in  the  analysis  of  infinite  arrays  was  done  by  Wheeler  and  Edelberg  and  Oliner, 

but  Stark10  gave  the  first  definite  exposition  of  the  complete  series.  Both  Wheeler11  and  Parad12 

13 

have  written  on  the  physical  interpretation  of  the  result,  and  Rhodes  has  given  an  elegant  der¬ 
ivation  of  a  related  expression  for  single  elements.  The  derivation  is  straightforward,  but  will 
only  be  sketched  here.  It  is  based  on  the  fact  that  from  an  expression  for  the  total  power  flow 
(both  real  and  reactive)  outward  from  a  typical  cell  of  an  infinite  array  at  the  array  face,  the 
driving  impedance  the  cell  presents  to  a  single -mode  transmission  line  can  be  calculated  and 
shown  to  be  related  to  the  element  pattern. 

To  illustrate,  let  us  consider  an  array  of  elements  on  a  rectangular  grid.  The  elements 

are  assumed  to  be  identical  metallic  conductors  of  infinite  conductivity  with  a  known  surface 

current  density  distribution  in  the  y-direction,  K  (x,  y).^  The  surface  current  density  is  identi- 

y  th 

cal  from  cell  to  cell,  except  that  the  relative  phase  of  the  current  at  the  center  of  the  mn  cell 
is  exp(-j(ma  +  n/3)].  Thus,  the  total  surface  current  density  distribution  can  be  expressed  as 
a  doubly  infinite  series: 


K(x,y)  =  ey  Y,  Z  Ky<x  _  mDx<  y  ~  nDy)  exp[-j(ma  +  n/3)] 
m  n 


(12) 


where  ey  is  the  y-directed  unit  vector.  Note  that  K(x,  y)  is  not  a  periodic  function;  however,  the 

product  K(x,  y)  exp[j(«x/D  +  /3y/D  )]  is  periodic  and  can  be  expanded  in  a  two-dimensional  Fourier 

x  y 

series.  The  result  is  that  the  normalized  current  density  can  be  written  as 


=  %  Z  E  Ky  <«.*>  exp[-j(kx  x  +  k  y)] 


me 


p  q 


pq 


where  we  have  defined 


k  = 
x  D 

P  x 

k  =  L+2 3£ 

yP  I5y 


and 


(13a) 

(13b) 


(13c) 


t  Only  th*  requirement  of  identical  element!  it  necessary;  the  other  restrictions  are  for  convenience  and  can  be 
removed. 
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The  definition  of  k?  is  included  for  later  reference,  and  kQ  is  the  free -space  propagation  con- 
2  2  M 

stant,  kQ  -  w  A'0«0.  By  the  usual  procedure,  the  Fourier  coefficients  are  found  to  be 


i  rV2  rV2  Kv(x*y) 

ypq  xDy  J-D^2  J-D^2  xinc 


exp  [j(k  x  +  k  y)]  dxdy 
P  yq 


If  we  now  apply  Maxwell's  equations  to  the  free -space  region  in  front  of  the  array,  we  obtain 
expressions  for  the  electric  and  magnetic  fields  in  space  in  terms  of  the  surface  current  density 
at  the  aperture.  The  electromagnetic  fields  in  space  must  be  periodic  in  the  same  sense  as  the 
surface  current  density  and,  in  fact,  must  be  an  infinite  set  of  propagating  or  attenuating  plane 
waves.  For  the  electric  and  magnetic  fields  tangential  to  the  aperture  plane,  we  have 


Ky  (a ,  0) 

H  (x,  y,  z)  =22  — ^ -  exp[-j(k  x  +  k  y  +  k  z)] 

p  q  P  yq  pq 


(15a) 


E  (x,y,z)  =  ^  2  2 


<ko2-kv)Kv  <"■*> 

J  rt  *  r\n 


exp  [— j(k  x  +  k  y  +  k  z))  (15b) 

XP  yp  zpq 


where  £  =  J P0/c0  is  the  intrinsic  impedance  of  free  space. 

It  should  be  noted  that  the  wave  equation  for  the  free-space  region  defines  an  infinite  spec¬ 
trum  of  plane  waves.  Each  plane  wave  is  associated  with  a  "propagation  vector"  with  compo¬ 
nents  k  ,  k  ,  and  k  .  In  general,  these  components  can  be  complex,  but  for  this  case,  k  and 
x  y  z  x 

ky  are  constrained  to  be  real,  while  k^  is  positive  real  or  negative  imaginary  such  that  the  wave 
equation  will  be  satisfied: 


k2  =  k2-k2-k2 
z  o  Kx  Ky 


2  _ 
ui  |i  t 
o  o 


Physically,  for  k  2  +  ky  <  k^ ,  we  can  identify  the  k's  with  the  angles  of  Fig.  3: 

k  =  k  sin  9  cos  <p 
x  o 

k  =  k  sin  9  sin  <t> 
y  O 

k  =  k  cos  9  (16) 

z  o 

In  this  notation,  for  example,  the  far  field  of  a  typical  element  with  normaliz  ed  current  density 
[Ky(X*y)/linc>  ey  is 


r  D  /  2  fD/2  K  (x,  y) 

\  x  i  y  _y__ —  exp  [j(k  x  +  kv))  dxdy 

v  -  D  / 2  J-D ^2  Xinc  x  r 


where  I.  is  the  total  incident  current  in  a  cell.  The  usual  pattern  in  and  9  coordinates  is 
inc 

obtained  by  substitution  of  Eqs.(16)  into  this  result.  Note  that  the  Fourier  components  of  the 

2  2  2 

current  density  are  related  to  the  far-field  pattern  of  an  element  when  k  +  k  ^  k  ,  Although 

y  _  o 

222  p  q 

we  cannot  associate  the  plane  waves  that  have  kx  +  ky  >  kQ  with  the  element  radiation  pattern, 

they  still  have  significance  because  they  correspond  to  stored  energy  in  the  aperture. 

It  is  also  worth  noting  that  the  discrete,  periodic  boundary  condition  defined  by  the  aperture 

of  the  array  causes  the  radiated  power  to  be  confined  to  a  discrete  set  of  directions  in  k  k  -space. 

x  y 
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In  the  above  notation,  this  means  that  the  variables  k  and  k  take  on  only  the  discrete  set  of 

x  y 

values  k  and  k  for  every  beam-pointing  angle  specified  by  the  element  phasings,  a  and  0. 

xp  yq 

This  contrasts  with  the  situation  for  finite,  planar  apertures  where  k  and  k  are  continuous 

x  y 

variables. 

If  the  normal  component  of  the  Poynting  vector,  -E  H*,  is  integrated  over  one  cell  in  the 

y 

aperture  plane,  we  obtain  an  expression  for  the  complex  radiated  power  per  element.  The  ele¬ 
ment  driving  impedance  is  related  to  the  radiated  power  and  is  given  by 

dd  00  *  (k  2  -  k  2  ) 

ZD(a>/3)  *  *  -fa*  2  2  ~k— ^  |Kv  •  (17) 

°  P=~  Zpq  Pq 

Observe  that  the  numerator  in  the  p  =  q  =  0  term  is  just  the  power  pattern  of  an  element  of  the 
array  when  (a/Dx)2  +  (/?/Dy)2  <  k  2. 


1.  Physical  Interpretation  of  the  Series 


To  grasp  the  significance  of  the  "directions"  specified  by  the  factors  k  ,  k  ,  and  k  , 

xp  yq  pq 

note  that  we  can  write  the  array  factor  of  a  uniformly  illuminated  array  with  phase  constants  a 
and  0  as 

A(kx,  ky)  *  £  Y,  exp  [j(ma  +  n/3)]  exp[-j(kxmDx  +  kynDy)] 
m  n 


As  it  stands,  A  is  a  function  of  two  variables  kx  and  ky  that  can  range  from  -*®  to  +<*>.  It  will 
have  maxima  at  points  "in  k-space"  where 

(kxDx  -  a)  =  pin 

(kyDy  “  p)  =  q2x 

or,  referring  to  Eq.  (17),  at  all  values  of  k  and  k  satisfying 

x  y 

kx  =  kx  ,  p  =  0,  ±1,  ±2, . . . 

P 

k  =  k  ,  q  =  0,  *1,  ±2,  .  .  . 
y  yP 

That  is,  the  infinite  array  generates  a  "nest"  of  delta  functions  in  k-space  on  a  regular  grid  as 
indicated  in  one  dimension  in  Fig.  6(a).  The  delta  functions  inside 


k2+k2=k2 
x  y  o 


correspond  to  grating  lobes  in  visible  space,  the  remainder  to  lobes  in  "invisible  space." 

Consequently,  we  have  a  physical  interpretation  of  each  term  of  Eq.  (17):  the  pq  term  cor¬ 
responds  in  amplitude  to  the  element  "power  pattern"  {Fig.  6(b))  of  a  typical  element  in  the  "di¬ 
rection"  k  ,  k  divided  by  the  z-direction  "propagation"  constant  (Fig.  6(c)).  Fork  and  k 

P  yQ  Xp  yq 

in  real  space,  the  z-direction  propagation  constant  is  real  and  that  term  of  the  series  contributes 

only  to  the  real  part  of  the  impedance.  If  k  and  k  are  large  enough  to  make  k  imaginary, 

p  yq  zpq 

that  entire  term  will  be  imaginary  and  will  contribute  only  to  the  reactance. 

An  important  observation  about  the  symmetry  of  the  driving  impedance  variation  with 
scan  follows  from  the  fact  that  substitution  of  -a  and  -0  into  Eqs.  (13)  for  k  and  k  changes 

xp  yq 
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Fig.  6.  Factors  of  grating-lobe  series  (one  dimensional):  (a)  scanning 
array  of  delta  functions,  (b)  element  pattern,  (c)  1/k  factor,  (d)  re¬ 
sulting  contributions  to  series  for  Z^(8). 

etc.  Since  we  sum  over  all  p  and  q,  the  signs  of  p  and  q  are  unimportant. 


Further,  since  K  (a ,  0)  is  multiplied  by  its  complex  conjugate  in  Eq.  (17),  the  sign  ZD(or,/3) 

ypq 

is  invariant  to  changes  in  sign  of  k  and  k  ;  hence, 

P  yq 

ZD(<*,0)  =  Z  D(— =  Z  D(a,  — /?)  =  ZD(-a,-0)  (18) 

for  an  element  irf  an  infinite  array.  Thus,  in  the  following  examples,  we  will  only  examine  the 
behavior  of  ZD(a,/J)  for  positive  a  and  ft. 

2,  Some  Qualitative  Inferences  from  the  Series 

The  effect  of  scanning  the  beam  on  the  driving  impedance  of  the  element  thus  depends  on  the 
change  of  the  sum  of  the  contributions  of  the  individual  grating  lobes  as  the  nest  of  grating  lobes 
is  scanned  in  k-space.  For  example,  the  effect  of  a  grating  lobe  crossing  the  boundary  between 
real  and  imaginary  space  depends  markedly  on  the  behavior  of  the  element  pattern  at  that  bound¬ 
ary,  since  l/kz  =  °° .  Unless  the  element  pattern  has  a  null  in  that  direction,  a  discontinuity  in 
impedance  occurs  as  the  grating  lobe  crosses  over  this  boundary,  and  abruptly  changes  its  con¬ 
tribution  from  reactance  to  resistance  or  vice  versa. 

The  limiting  cases  of  the  effect  of  element  spacing  on  the  element  impedance  variation  can 
be  easily  assessed  from  the  grating-lobe  series.  In  the  limiting  case,  as  the  spacing  approaches 
zero,  the  grating  lobes  recede  toward  infinity.  As  the  elements  shrink  to  accommodate  the  spac¬ 
ing,  the  current  density  becomes  a  constant  for  all  cells,  and  the  element  power  pattern,  the 

2  /  2 

numerator  in  the  summation  of  Eq.  (17),  is  just  the  "obliquity  factor,"  1  -  ky/kQ  .  The  pq  term 

of  the  scries  thus  has  the  following  asymptotic  behavior  for  large  k  and  k  (small  D  and  D  ) 

P  yq  Y 


k2-k2 

TTTT 

o  z 

pq 


-jk  along  the  k  axis 

yq  y 

along  the  kx  axis 
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and  consequently  the  reactance  becomes  very  large  due  to  the  impulses  at  large  values  of  k  . 
However,  the  variation  in  impedance  due  to  the  lobes  other  than  p  =  q  =  0  tends  to  zero,  since 
the  amount  of  phase  differential  P/®y  between  elements  required  to  steer  a  beam  to  a  given 
angle  is  constant  and,  therefore. 


lim  k 
D.  -*0  yc 


2ir 

Dy 


which  is  independent  of  element  phasing.  Only  the  impedance  contribution  of  the  p  =  q  =  0  term 
varies,  and  the  fractional  variation  in  is  seen  from  Eq.  (17)  to  be 

RP(Q)  ^  1  -sin2  6  sin2 

r  '  Rd(0,  0)  cos  0 

For  E-plane  scan,  <f>  =  ir/2  and 
r „  =  cos  0 


while  for  H -plane  scan,  <f>  =  0  and 

1 

rH  =  cos  0 

in  agreement  with  Eqs.  (8)  and  (11)  for  the  resistive  sheet  and  the  trend  indicated  by  the  examples 
of  Fig.  5,  if  account  is  taken  of  the  fact  that  Eqs.  (8)  and  (11)  are  for  the  impedance  looking  into 
the  array  from  space.  The  above  results  are  for  the  impedance  looking  out.  The  normalized 
impedance  is  thus  inverted. 

For  widely  separated  elements  of  small  active  area,  the  grating-lobe  series  indicates  ex¬ 
actly  what  intuition  suggests:  as  the  element  separation  increases,  the  impulses  move  very  close 
together  so  that  many  are  in  visible  space,  and  scanning  the  array  merely  moves  some  out  of 
visible  space  and  others  into  it,  with  a  net  result  of  little  impedance  variation  with  scan.  Of 
course,  as  a  consequence,  there  are  also  numerous  lobes  in  the  actual  spatial  pattern  in  visible 
space,  as  indicated  in  Fig.  7,  with  resulting  loss  of  sensitivity  and  potential  directional  ambiguity. 


Fig.  7.  Effect  of  element  size  on  array  patterns: 
(a)  element  pattern,  (b)  array  factor  for  element 
spacing  D,  (c)  far-fie Id  pattern. 
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Of  more  practical  interest  for  elements  with  widely  separated  centers  is  the  case  in  which 
the  elements  are  themselves  directive.  One  way  to  achieve  this  is  to  have  the  current  distribu¬ 
tion  occupy  most  of  the  space  between  the  elements,  such  as  using  parabolic  reflectors.  The 
pattern  result  is  indicated  in  Fig.  7:  using  more  of  the  allotted  space  improves  the  suppression 
of  the  grating  lobes.  In  the  limit,  as  we  fill  all  allotted  space  with  a  uniform  current  distribu¬ 
tion,  the  pattern  grating  lobes  are  completely  suppressed  when  the  beam  is  pointed  at  broadside 
However,  they  are  only  partially  suppressed  as  ihe  beam  is  scanned.  The  same  thing  is  taking 
place  in  the  grating-lobe  series;  the  wide  center-to-center  element  spacing  places  many  of  the 
nest  of  impulses  in  visible  space,  but  these  are  suppressed  to  a  degree  which  depends  on  the 
extent  to  which  the  aperture  is  filled. 

3.  Additions  to  the  Basic  Series 


Several  interesting  and  useful  additions  can  be  made  in  the  basic  grating-lobe- series  for¬ 
mulation.  First,  as  Stark  has  pointed  out,  the  modification  of  Eq.  (17)  that  must  be  made  to  ac¬ 
commodate  elements  above  an  infinite,  perfectly  conducting  ground  plane  consists  of  the  addition 
of  a  factor  dependent  on  the  element-to-ground-plane  spacings: 


Z  =  | 


D  D 
x  y 


Z 

p- 


£ 


<k2-k2> 
o  y  ' 

‘'q 

pp 


I  K  (<*,0)1  { 1  —  exp  [— j2k  s)} 


PP 


PP 


(19) 


This  additional  factor  is  of  no  small  importance,  since  it  goes  to  zero  whenever  k  =  0  in  the 

zpq 

denominator.  The  effect  on  impedance  variation  is  pronounced,  as  we  shall  see. 

The  grating- lobe  series  is  also  adaptable  to  any  regular  grid  of  elements  (e.g.,  triangular). 
It  is  only  necessary  to  set  the  relative  spacings  and  placement  of  the  nest  of  impulses  to  cor¬ 
respond  to  the  grid  shape  and  spacing. 

The  series  can  also  be  used  to  infer  the  approximate  average  (with  respect  to  position  in  the 
array)  impedance  behavior  of  the  elements  of  small  arrays  of  identical  elements.  The  infinite 
array  assumption  in  the  derivation  is  responsible  for  two  aspects  of  the  results  of  Eqs.  (17)  and 
(19): 


(a)  Justification  for  assuming  that  all  elements  have  the  same  impedance 
and  known  relative  amplitude  and  phase  weighting  (see  Sec.II-A).  There¬ 
fore,  the  impedance  of  any  particular  element  is  simply  related  to  the 
variation  in  radiated  and  stored  power  of  the  entire  array. 

(b)  The  sampling  functions  (the  nest  of  delta  functions  of  Fig.  6)  have  infinites¬ 
imal  width. 


Thus,  to  apply  the  grating-lobe  series  to  small  arrays,  we  must  be  willing  to  settle  for  an  approx¬ 
imate  average  behavior  of  the  elements  (approximate,  since  we  still  assume  in  the  derivation 
that  the  relative  amplitudes  and  phases  of  the  elements  are  known,  which  actually  depend  on  the 
relative  impedances).  If  we  accept  this  shortcoming,  the  effect  of  the  small  array  is  to  replace 
the  delta  functions  with  sampling  functions  of  a  width  equal  to  the  array  beamwidth.  The  net  ef¬ 
fect  is  that  a  reasonable  qualitative  assessment  of  the  small  array  impedance  variation  can  be 
obtained  from  that  of  the  infinite  array  by  averaging  the  infinite  array  data  over  a  beamwidth  of 
the  small  array. 

Finally,  the  series  can  be  extended  to  include  arrays  of  crossed-dipole  pairs  with  or  with¬ 
out  ground  planes.  Since  the  orthogonal  dipoles  are  radiatively  coupled  to  each  other  at  nearly 
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all  scan  angles,  the  analysis  of  this  type  of  array  is  more  complicated.  That  is,  the  radiative 
coupling  will  influence  the  relative  currents  on  the  dipoles  of  every  pair  with  the  result  that  the 
type  of  feed  network  behind  every  dipole  pair  must  be  included  in  the  analysis.  The  problem  can 
still  be  formulated  in  rather  general  terms  for  arbitrary  three-  or  four-port  feed  networks  behind 
every  dipoie  pair;  however,  we  must  completely  specify  the  parameters  (scattering  matrix,  for 
example)  of  the  feed  network  before  valid  numerical  results  can  be  obtained.  In  the  analysis  of 
this  type  of  array,  we  must  consider  the  polarization  properties  of  the  transmitted  waves,  the 
impedance  variation  at  the  input  port  (or  ports)  of  the  feed  network  for  a  dipole  pair,  and  the 
power  dissipated  in  the  feed  network  or  its  termination. 

E.  Dependence  of  Element  Impedance  on  Element  and  Array  Design  Parameters 

In  this  section,  we  present  results  of  the  dependence  of  the  element  driving  impedance 
Z p(a,0)  on  the  pointing  direction  of  the  main  beam  of  the  array.  We  first  examine  in  a  planar 
array  context  the  extent  to  which  the  impedance  variation  can  be  expected  to  depend  on  element 
type  and  spacing  and  other  design  details.  We  then  examine  the  difference  in  behavior  in  linear 
arrays,  including  a  brief  examination  of  the  frequency  dependence  of  the  driving  impedance. 

1.  Planar  Array  Impedance  Behavior 

To  make  the  data  more  directly  useful  for  design,  we  will  discuss  ZD  for  planar  arrays  as 
a  function  of  the  scan  angles  in  each  of  three  planes  through  array  broadside.  The  symbols  0^,, 
Gjj,  and  ©D  will  be  used  to  define  the  angle  of  the  beam  from  broadside  in  the  E-plane,  H-plane, 
and  "diagonal  (D)"  plane  (45°  from  both  E  and  H),  respectively.  If  the  elements  are  polarized  in 
the  y-direction,  we  can  make  the  correspondence  between  a,  p,  and  0^,,  0^,  and  0^  by 

—  2ttD 

a  = — ^ — —  sin©H  ,  H-plane  scan 
-ZvD 

P  = - ^ ^  sin  ©E  ,  E-plane  scan 

y—z  7F  2jr 

— y  +  -^-5  =  -r  sin©.-.  ,  D-plane  scan 

x  y 

The  data  presented  were  computed  from  Eqs.  (17)  and  (19),  as  appropriate,  using  a  digital 
computer  to  sum  several  hundred  terms.  To  provide  a  crude  assessment  of  the  dependence  of 
impedance  variation  on  element  type,  three  significantly  different  current  distributions  within 
a  cell  were  examined: 

(1)  An  idealized  short  dipole;  infinitesimal  diameter,  effective  length  0.1A, 

(2)  An  idealized  half-wavelength  dipole;  infinitesimal  diameter,  K  (y)  = 

Ko  cos  try/ A,  y 

(3)  A  uniform  current  sheet  filling  the  entire  D  by  D  cell. 

x  y 

Models  (1)  and  (2)  represent  idealized  models  of  both  dipoles  and  their  complement,  slots  in  an 

infinite  ground  plane.  Model  (3)  is  a  crude  but  tractable  model  of  the  dual  of  an  open-waveguide 

radiator.  It  has  recently  been  demonstrated  experimentally  that  model  (2)  accurately  predicts 

14 

the  behavior  of  real  dipole  arrays  even  for  lengths  not  particularly  close  to  a  half  wavelength. 
Some  experimental  data  on  open  waveguide  array  behavior  are  presented  later. 
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We  firiSrt,  compare  the  impedance  behavior  of  these  three  models  in  a  common  configura¬ 
tion  -  a  planar  afz^ty  with  elements  on  a  0.55X  square  grid  (maximum  scan  without  grating  lobes: 
55°  in  E-  and  H -planes).  We  then  examine  the  modification  of  the  dipole  behavior  when  a  ground 
plane  is  used.  The  effect  ^/  element  spacing  and  grid  shape  on  impedance  variation  is  assessed 
for  dipoles  above  ground,  an'd-{he  interplay  between  element  spacing  and  ground  plane  height  is 
examined  as  an  indication  of  the  intent  to  which  the  designer  can  optimize  performance.  The 
behavior  of  linear  arrays  is  briefly  etxpmined  and  compared  with  planar  array  behavior.  Finally, 
we  present  an  example  of  the  frequency'dependence  of  the  driving  impedance. 

For  uniform  presentation  of  the  data  in  a  useful  form,  ZD  vs  angle  is  displayed  on  Smith 
charts  (see  Fig.  8).  We  present  normalized  Z^s  i.e.,  for  ZD(0)  =  RD(0)  +  jXD(0), 


Z  (0)  -jXD(0) 
ZD(G)norm  =  RD(0) 


(20) 


where  0  refers  to  the  value  of  ZD  when  the  beam  is  pointed  at  broadside.  While  not  the  only 
choice  of  interest,  there  are  three  reasons  for  this  format: 

(a)  It  displays  the  impedance  variation  as  it  would  occur  if  the  array  were 
matched  at  broadside  —  a  condition  that  maximizes  gain  (but  does  not 
minimize  impedance  variation). 

(b)  It  makes  the  results  easier  to  correlate  with  experimental  measure¬ 
ments;  the  effects  of  any  actual  matching  networks  in  the  antenna  are 
accounted  for  by  normalizing,  except  for  an  arbitrary  rotation  of  the 
Smith  chart. 

(c)  In  some  cases,  the  grating-lobe  series  does  not  converge  to  a  unique 
ZD,  but  does  converge  to  a  unique  value  of  [ZD(8)  -  jX^OJl/R^O). 


To  aid  the  reader,  for  model  (2),  which  is  used  for  comparison  of  different  array  configurations, 
we  will  include  ZD<0)  and  also  the  impedance  of  a  single  isolated  element  (above  a  ground  plane 
if  appropriate),  ZR  (see  Fig.  9).t 

2.  Comparison  of  Element  Types 


Figure  8  shows  the  behavior  of  a  typical  element  of  an  array  of  short  dipoles  (SD)  [for  short 
slots,  interchange  E  for  II  and  impedance  for  admittance].  The  most  striking  fact  is  the  supe¬ 
riority  of  the  behavior  in  the  E-plane  over  the  other  planes.  Unfortunately,  the  E-plane  is  atyp¬ 
ical,  since  the  factor  (k^  -  k^  )  has  a  zero  in  the  E -plane  at  the  same  angle  as  the  zero  of  k  in 

yq 


the  denominator  of  Eq.  (17). 
chart  whenever 


For  all  other  scans,  the  impedance  goes  to  the  edge  of  the  Smith 


(a)  A  grating  lobe  comes  into  visible  space,  which  occurs  at  55"  in  the 
H -plane,  or 

(b)  The  main  lobe  is  scanned  out  of  visible  space,  which  occurs  at  90*  in 
the  D-plane. 

The  E-  and  H  -plane  curves  are  discontinued  at  65*.  This  angle  physically  corresponds  to  a 
grating  lobe  at  -65*  (symmetrically  disposed  relative  to  the  main  beam),  so  that  further  scanning 
is  equivalent  to  scanning  the  main  beam  in  from  65*;  the  curve  retraces  itself. 

If  we  adopt  as  an  arbitrary  criterion  a  tolerable  VSWR  of  3:1,  we  can  scan  to  about  47*  in 
the  H-plane,  about  79*  in  the  D-plane,  and  anywhere  in  the  E-plane;  for  a  2:1  VSWR,  40*  in  H, 
76*  in  D,  and  anywhere  in  E. 


t  For  wide  icon  angles,  a  match  to  Zg  very  nearly  minimizes  the  miimateh  with  scan,  but  at  a  loss  in  broodfida 
pain  with  respect  to  a  match  to  Z-.(0). 


Fig*  10.  Normalized  impedance  vs  scan  angle  in 
E-,  H-,  and  D -planes  for  a  typical  element  of  an 
array  of  uniformly  illuminated  cells;  D  =  D  = 

o.55x.  x  y 


Fig.  12.  Normalized  impedance  vs  scan  angle  in 
E-,  H-,  and  D-planes  for  a  typical  element  of  an 
array  of  short  dipoles  (l  *  0.  I  X)  mounted  0. 25  X 

above  ground  plane;  D  *  D  ■  0. 55  X. 

x  y 


Fig.  II.  Comparison  of  measured  impedance  variation 
of  a  central  element  ino  10  X  10  open-waveguide  array 
with  that  computed  for  the  same  element  in  an  identical 
array  of  X/2  dipoles  without  ground  plane;  D  =  D  = 

o.6x.  x  y 


Fig.  13.  Normalized  Impedance  vs  scan  angle  In 

E-,  H-,  and  0-planes  for  a  typical  element  of  an 

array  of  V?  dipoles  mounted  0. 25  X  above  ground 

plane;  D  *  D  *0.55X. 
x  y 


For  a  half-wavelength  dipole  (HD),  a  s  indicated  in  Fig.  9,  the  behavior  is  qualitatively  quite 
similar,  suggesting  that  the  dipole  length  is  of  secondary  importance.  For  a  3:1  VSWR,  we  can 
scan  to  about  45*  in  H,  79*  in  D,  and  anywhere  in  E;  for  a  2:1  VSWR,  about  40#  in  H,  77*  in 
D,  and  50*  in  E. 

The  "open-waveguide"  (OWGD)  dual  model  shows  a  still  greater  E-plane  variation  than  either 
dipole,  but  is  again  qualitatively  similar  (an  observation  made  more  obvious  by  comparing  the 
OWGD  dual  model  results  of  Fig.  10  with  Fig.  9  rotated  about  45°  ccw).  For  a  3:1  VSWR,  we  can 
scan  about  51*  in  H,  about  75°  in  D,  and  about  35°  in  E;  for  a  2:1  VSWR,  about  42°  in  H,  45° 
in  D,  and  about  28*  in  E. 

These  results  are  summarized  in  Table  I. 


TABLE  1 

IMPEDANCE  BEHAVIOR  OF  ELEMENTS 
WITHOUT  GROUND  PLANES 

Dement  Type 

VSWR 

Scon  Plane 

SD 

HD 

OWGD 

E 

35 

3:1 

H 

47 

45 

51 

D 

79 

79 

75 

E 

— 

50 

28 

2:1 

H 

40 

40 

42 

D 

76 

77 

45 

Since  the  OWGD  dual  model  is  crude,  its  apparent  inferiority  should  not  be  taken  too  seri¬ 
ously.  To  reinforce  this  point  of  view.  Fig.  11  shows  a  comparison  of  experimentally  determined 
Zj-j  for  a  finite  (10  x  10)  array  of  open -waveguide  radiators  with  computations  for  a  similar 
(10  x  10)  array  of  half- wavelength  dipoles  without  a  ground  plane.  In  this  case,  the  impedance 
variation  of  the  element  in  the  waveguide  array  is  seen  to  be  slightly  less  than  that  for  the  same 
element  in  the  dipole  array  for  scan  in  the  waveguide  E-plane  (the  H-plane  for  dipoles).  For 
scan  in  the  waveguide  H-  and  D-planes,  however,  the  results  for  the  two  arrays  are  within  the 
limits  set  by  experimental  tolerances.  Similar  results  were  obtained  for  other  elements  in  this 
particular  array  and  for  waveguide  and  dipole  arrays  with  different  inter-element  spacings 

(D  -  0.6 X,  D  =  0.3X).  The  obvious  conclusion  is  that  the  impedance  behavior  of  an  element  in 
x  y 

a  waveguide  array  can  be  predicted  reasonably  well  by  computations  on  a  similar  array  of  dipoles. 

The  addition  of  a  ground  plane  under  the  dipoles  has  a  marked  effect  on  the  impedance  be¬ 
havior  as  shown  in  Figs.  12  and  13  for  the  short  and  half-wavelength  dipoles.  Physically,  the 
addition  of  the  ground  plane  below  an  electric  dipole  prevents  the  dipole  from  radiating  along 
the  ground  plane;  hence,  in  the  language  of  the  grating -lobe  series  of  Eq.  (19),  the  ground  plane 
factor,  1  -  exp(-j2k  s],  approaches  zero  as  k  —  0,  so  that  the  impedance  is  continuous  when  a 
grating  lobe  becomes  visible.  Also,  if  an  element  can  radiate  in  the  direction  of  the  other  ele¬ 
ments,  a  grating  lobe  directed  toward  the  other  elements  will  cause  them  to  absorb  as  much 
power  as  they  radiate,  bringing  about  a  unity  reflection  coefficient. 
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The  diagonal  plane  scan  still  causes  the  reflection  coefficient  to  reach  a  magnitude  of  unity, 
since  the  beam  can  be  scanned  out  of  visible  space  without  a  grating  lobe  coming  into  visible 
space,  thus  suppressing  all  radiation  (for  an  infinite  array). 

For  the  elements  mounted  a  quarter -wavelength  above  ground,  the  allowable  scamj  for  3:1 
and  2:1  VSWR's  are  as  shown  in  Table  II. 


TABLE  II 

IMPEDANCE  BEHAVIOR  OF  ELEMENTS 
WITH  GROUND  PLANES 

Element  Type 

VSWR 

Scan  Plane 

SD 

HD 

E 

55 

50 

3:1 

H 

50 

50 

D 

62 

62 

E 

45 

40 

2:1 

H 

40 

40 

D 

52 

50 

The  differences  between  the  short-  and  half-wavelength  dipole  array  results  are  not  significant, 
particularly  in  view  of  the  ability  to  alter  the  impedance  behavior  somewhat  by  varying  s,  as  we 
shall  see. 

In  summary,  the  most  striking  difference  in  impedance  behavior  of  the  elements  investigated 
is  brought  about  by  suppressing  the  ability  of  an  element  to  radiate  in  the  direction  oi  other  ele¬ 
ments.  For  dipoles  and  some  similar  elements,  this  can  be  easily  accomplished  by  placing  them 
above  a  ground  plane.  Note,  however,  that  slots  or  open  waveguides  in  a  ground  plane  do  radi¬ 
ate  along  the  ground  plane  and  hence  are  analogous  in  behavior  to  a  dipole  without  a  ground  plane 
(the  thin  slot  in  a  ground  plane  is  the  exact  dual  of  the  dipole  without  a  ground  plane). 

Even  this  difference  in  behavior  is  only  important  for  scan  angles  near  grating-lobe  forma¬ 
tion.  It  is  evident  from  Tables  I  and  II  that  if  the  usual  element  spacing  criterion  of  allowing  a 
grating  lobe  to  appear  marginally  is  used  to  design  the  array,  elements  such  as  dipoles  would 
give  vastly  superior  (but  still  poor)  VSWR  performance  at  extreme  scans  if  placed  above  a  ground 
plane.  If,  however,  the  element  spacing  is  chosen  to  constrain  the  maximum  VSWR  to  a  reason¬ 
able  value,  there  appears  to  be  little  difference  between  the  elements,  and  other  considerations 
(e.g.,  mechanical  convenience)  should  dictate  the  choice.  Note,  however,  that  such  a  criterion 
results  in  a  much  smaller  allowable  scan  volume  for  a  given  element  grid  spacing. 

To  illustrate  the  extent  to  which  detailed  design  parameters  may  affect  impedance  behavior, 
we  will  explore  the  interlocking  effects  of  two  parameters  of  an  array  of  dipoles  above  ground. 

We  will  also  examine  the  effect  of  the  height  of  the  dipole  above  ground,  s,  on  the  maximum 
VSWR  which  occurs  for  a  specified  scan  volume.  In  particular,  we  will  choose  a  40*  cone- 
shaped  scan  volume  and  examine  elements  on  a  square  grid  (suggested  by  the  symmetry  of  the 
scan  volume)  and  a  triangular  grid. 

A  plot  of  maximum  VSWR  as  a  function  of  ground  plane  spacing  is  shown  in  Fig.  14  for  E- 
and  H-plane  scans  of  two  square-grid  arrays.  From  this  figure,  it  is  apparent  that  the  maximum 
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Fig.  14.  VSWR  at  6max  vs  ground  plane  spacing  g 
for  large,  square-grid  arrays  of  )s/2  dipoles. 
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E-plane  mismatch  becomes  very  large  and  the  maximum  H-plane  mismatch  becomes  monoton- 
ically  smaller  as  the  ground  plane  spacing  is  decreased.  Thus,  for  every  inter-element  spac¬ 
ing,  there  is  an  optimum  ground  plane  spacing  which  equalizes  the  maximum  mismatch  in  the 
principal  planes  of  scan  for  a  specified  scan  volume. 

A  similar  set  of  computations  for  planar  arrays  with  the  elements  on  an  equilateral  trian¬ 
gular  grid  is  shown  in  Fig.  15  for  a  scan  of  40°.  Qualitatively,  the  results  are  quite  similar  to 
those  for  the  square -grid  planar  arrays.  In  particular,  for  every  inter-element  spacing,  there 
is  an  optimum  element -to -ground -plane  spacing  which  minimizes  the  maximum  mismatch  (as¬ 
suming  a  match  at  broadside)  incurred  for  a  conical  scan  volume. 

Figure  16  shows  a  plot  of  impedance  variation  with  scan  angle  for  a  triangular-grid  array 
mounted  at  the  optimum  height  above  the  ground  plane  (compare  with  Fig.  13  for  the  square-grid 
array). 

A  cursory  examination  of  the  data  for  square-grid  and  triangular-grid  arrays  with  approx¬ 
imately  equal  areas  per  element  indicates  that  the  triangular -grid  array  may  be  preferable  on 
the  basis  of  impedance  variation,  although  the  difference  is  slight.  Further,  the  comparison  is 
not  direct,  since  the  two  arrays  have  somewhat  different  spatial  coverages. 

For  all  the  linear-  and  planar-array  configurations,  it  is  seen  that  for  optimum  ground 
plane  spacing  for  a  given  volume  coverage,  an  array  with  elements  more  closely  spaced  than 
required  by  grating-lobe  considerations  will  exhibit  less  impedance  variation  than  one  which 
just  satisfies  the  scan  volume  requirement:  this  is,  of  course,  obtained  at  the  expense  of  a 
larger  number  of  elements  to  realize  the  same  gain  and  beamwidth. 

3.  Linear  Array  Impedance  Behavior 

For  a  planar  array  with  different  scan  requirements  in  the  x-  and  y-directions,  an  additional 
design  parameter  is  the  choice  of  element  alignment  with  the  x-  and  y-axes.  This  is  most  evi¬ 
dent  in  what  we  might  view  as  an  extremely  degenerate  case:  a  linear  (one-dimensional)  array. 
Furthermore,  linear  arrays  are  obviously  of  interest  in  their  own  right.  Consequently,  we  will 
briefly  examine  the  difference  in  impedance  behavior,  using  arrays  of  dipoles  above  ground  for 
illustration. 
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Fig.  17.  Normalized  impedance  vs  scan  angle 
for  typical  elements  in  large,  linear  arrays  of 
X/2  dipoles  mounted  0.25  X  above  ground  planes; 
D  =  0.6  X. 


The  relative  impedance  variation  with  scan  angle  is  shown  in  Fig.  17  for  a  linear  array  of 
dipoles  spaced  0.6  X  on  centers,  0.25  X  above  a  ground  plane  for  three  orientations  of  the  dipoles 
relative  to  the  direction  of  the  array.  Observe  that  nearly  pure  resistance  variation  with  scan 
angle  is  obtained  when  the  dipoles  are  mounted  collinearly,  while  for  parallel  dipoles,  the  varia¬ 
tion  is  mainly  reactive.  Also  note  that  the  mismatch  at  the  maximum  scan  angle  is  approximately 
the  same  for  all  three  dipole  orientations.  Unfortunately,  neither  of  these  results  applies  for  all 
ground  plane  and  inter-element  spacings;  that  is,  the  above  choice  of  array  parameters  gives 
unique  results  and  is  presented  merely  to  emphasize  the  different  impedance  behavior  to  be  antic¬ 
ipated  as  a  function  of  polarization  in  a  linear  array. 

In  Fig.  18,  array  performance  is  assessed  as  a  function  of  inter-element  and  ground  plane 
spacings  (assuming  a  match  at  broadside)  by  giving  the  VSWR  incurred  at  the  maximum  scan 
angle  as  a  function  of  the  spacings.  From  this  figure,  it  is  readily  evident  that  a  greater  im¬ 
pedance  variation  with  scan  angle  is  exhibited  by  arrays  of  collinear  dipoles  for  most  ground 
plane  spacings.  The  variation  of  maximum  VSWR  with  ground  plane  spacing  (up  to  s  =  0.25  X)  is 
seen  to  be  small  for  arrays  of  parallel  and  45* -echelon  dipoles. 

If  the  frequency  is  now  allowed  to  vary  over  a  20-percent  band  for  an  array  with  D  =  0.6  X 
and  s  =  0.25  X,  the  impedance  variation  with  scan  angle  at  band  center  and  at  the  band  edges  is  as 
shown  in  Fig.  19.  Results  are  given  for  collinear  dipole  and  parallel  dipole  arrays.  The  fact 
that  the  impedance  variation  at  the  band  edges  is  not  substantially  greater  than  at  band  center 
indicates  that  some  other  component  (or  components)  in  the  feed  network  will  probably  limit  the 
system  bandwidth  more  than  the  antenna.  (It  is  also  interesting  to  observe  that  the  Smith  chart 
impedance  plots  for  the  parallel  dipole  and  collinear  dipole  arrays  remain  approximately  orthog¬ 
onal  as  the  frequency  is  varied.) 

ra.  COUPLING  IN  REGULAR  AND  IRREGULAR  FINITE  ARRAYS 

Intuitively,  we  would  expect  that  the  infinite  array  result  discussed  in  the  preceding  section 
would  accurately  describe  the  impedance  behavior  to  be  expected  for  elements  imbedded  in  even 
an  irregular  finite  array,  as  long  as  the  elements  in  the  neighborhood  of  an  element  in  question 
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provide  a  uniform  environment;  i.e.,  the  element  is  not  near  the  edge  of  the  array  and,  if  the 
array  is  irregular  in  spacing,  the  departure  from  regularity  is  small  over  a  region  of  several 
elements.  If  this  statement  is  to  be  useful,  we  must,  of  course,  have  quantitative  data  on  the 
degree  to  which  elements  in  nonuniform  surroundings  depart  from  the  infinite  array  prediction. 

We  shall  see  that  this  departure  depends  on  the  element  details,  the  location  of  the  other  ele¬ 
ments,  and  the  feed  network  that  is  used  to  excite  the  elements. 

Furthermore,  we  are  now  in  a  position  to  discuss  the  effects  of  coupling  on  array  patterns 
(the  patterns  of  infinite  arrays  are  singularly  uninteresting).  Here  again  the  nature  of  the  ele¬ 
ment  feed  is  important.  For  independently  fed  elements  (separate  excitation  or  an  isolating 
feed  network  such  as  a  well-matched  corporate  feed  using  four-port  junctions),  we  can  arrive 
at  some  very  general  results;  for  feed  networks  which  allow  element  interaction  through  the 
feed,  we  can  only  outline  a  method  of  analysis  and  present  some  representative  results  for  a 
few  specific  cases. 

A.  Methods  of  Analysis 

Since  the  grating-lobe  series  can  only  supply  an  average  (in  the  sense  explained  in  Sec.  II-D-3) 
element  impedance  behavior  in  a  finite  array,  the  first  problem  is  to  characterize  the  array  (in¬ 
cluding  mutual  coupling)  in  a  form  suitable  for  mathematical  analysis.  This  can  be  done  by  taking 
advantage  of  the  well-known  fact  that  the  terminal  currents  and  voltages  of  the  elements  in  an  ar¬ 
ray  can  be  represented  by  circuit  equations  as  Eq.  (1),  which  is  repeated  here  for  reference: 

V™=  ^m».Pq'pq  ' 

Since  there  are  N  antennas,  there  will  be  N  equations  of  the  form  of  Eq.  (1).  This  is  analogous 
to  the  case  of  an  N  terminal-pair  network.  The  set  of  equations  can  be  written  in  the  matrix 
form 

V)  =  [Z]  I]  (21) 

where  V]  and  I]  are  N  element  column  matrices,  and  [Z]  is  an  N  x  N  impedance  matrix. 

We  could  equally  well  relate  I  to  V  by  an  admittance  matrix  (inverse  of  [Z])  or  relate 

inn  PQ 

incident  and  reflected  waves  by  a  scattering  matrix.  Of  these,  the  scattering  matrix  is  the 
easiest  to  measure  in  practice  because  it  involves  only  a  measurement  of  the  power  (magnitude 
and  phase)  coupled  to  antenna  n  when  antenna  m  is  driven  with  unit  power.  Measurements  of 
impedance  or  admittance  coefficients  are  generally  much  more  difficult  because  they  involve 
either  open-circuit  voltages  or  short-circuit  currents. 

None  of  the  matrices  is  easy  to  compute  theoretically  for  most  types  of  array  elements,  be¬ 
cause  electromagnetic  scattering  (including  multiple  scattering),  as  well  as  induction  and  static 
field  coupling,  must  be  included  in  the  computation  of  a  single  coefficient.  However,  dipoles  of 
any  length  with  infinitesimal  diameters  can  be  handled  analytically  by  the  induced  emf  method 
for  computing  mutual  impedances*5  The  admittance  and  scattering  matrices  can  then  be  com¬ 
puted  from  the  impedance  matrix.  Figure  20  shows  the  behavior  of  the  mutual  impedances  be¬ 
tween  k/2  dipoles  0.25  k  above  ground  as  a  function  of  relative  position.  Further  data  for  dipoles 
and  slots  are  given  by  Kraus.  15 

The  second  important  part  of  the  analysis  of  finite  arrays  is  then  the  characterization  of  the 
network  feeding  the  antennas.  These  feed  systems  can  take  any  one  of  several  forms,  e.g.,  the 
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Fig.  20.  Mutual  impedance  between  dipoles  mounted  0.25  X  above 
ground  plane  as  a  function  of  relative  position. 
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Fig.  21.  Basic  serial  RF  array  feed  with  phase 
shifters  In  the  branch  lines. 
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F;g.  22.  Generator  equivalent  circuit 
for  arrays  with  independent  feeds. 
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serial  feed  of  Fig.  21.  Basically,  they  can  be  divided  into  two  categories:  independent  (isolating) 
feeds  and  nonisolating  feeds.  An  independent  feed  is  one  in  which  the  feed  to  an  antenna  in  the 
array  is  not  influenced  by  the  behavior  of  the  remainder  of  the  feed  system  except  through  radia¬ 
tion  coupling  between  antennas.  Examples  of  independent  feeds  are  those  constructed  with  direc¬ 
tional  couplers  and  those  using  circulators  or  isolators  behind  each  antenna.  A  nonisolating  feed, 
on  the  other  hand,  is  one  which  does  not  meet  the  above  criteria,  that  is,  a  feed  network  in  which 
there  is  internal  cross -coupling  within  the  structure. 

The  independent  source  problem  can  most  easily  be  analyzed  directly  from  the  impedance 
matrix  by  representing  the  sources  as  N  independent  generators  with  open-circuit  voltages  en 

and  impedances  Z  ,  as  indicated  by  the  equivalent  circuit  of  Fig.  22.  The  open-circuit  voltages 
gn 

are  assumed  to  have  independently  controllable  amplitudes  and  phases.  The  terminal  currents 
and  voltages  on  the  antennas  are  then  related  to  the  source  parameters  by 


V 

n 


n 


n  =  i,  2.  . . .  N 


or,  in  matrix  form, 

V)  =  e] -[Zg]  •  I]  (22) 

where  [Zg]  is  an  N  x  diagonal  matrix  with  the  diagonal  terms  given  by  the  N  generator  imped¬ 
ances  of  the  sources.  When  Eq.  (22)  is  substituted  in  Eq.  (1),  we  have 

e)  =  {[Z]  +  (Zg]>  •  I]  (23) 

Since  the  antenna  currents  are  unknown,  the  impedance  matrix,  modified  by  the  addition  of  the 
generator  impedances  to  the  antenna  self-impedances,  must  be  inverted  to  obtain  an  admittance 
matrix  ( Y]  relating  the  antenna  currents  to  the  generator  voltages 

I]  =  [Y1  e]  .  (24) 

Once  the  currents  have  been  computed  for  a  particular  set  of  drive  voltages  (a  particular 
beam-pointing  angle),  it  is  a  simple  matter  to  compute  the  element  driving  impedance  for  any 
one  of  the  antennas  from 


©)  =  T51  -Z 

m  *m 


(25) 


It  should  be  recalled  that  for  infinite  arrays,  we  can  rewrite  Eq.  (5)  in  terms  of  polar  angles 
(see  Fig.  3)  as 


Z  z- 

x  exp  {-jk[D  (m  —  p)  sin  0  cos  $  +  D  (n  -  q)  sin  0  sin  $)}  (26) 

x  y 


V'0,  =  Za+  -  “mn.pq 
p,  q?fem,n 


where  we  generalize  to  a  particular  location  x  =  mD^,  y  =  nD^  for  the  element  i.i  question.  From 
this  equation,  we  have  an  alternate  method  to  the  grating-lobe  series  for  computing  the  active  im¬ 
pedance  of  an  element  in  a  very  large,  regular,  uniformly  illuminated  array  without  resorting  to 
a  matrix  inversion. 

This  approximation  has  been  found  to  be  extremely  good  for  the  center  element  in  a  large 
dipole  array  above  a  ground  plane  (the  approximation  is  not  as  good  when  the  ground  plane  is 
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absent).  In  fact,  the  active  impedance  of  the  center  element  in  a  65  (collinear)  x  149  (parallel) 
array  (with  a  ground  plane)  agrees  within  0.1  percent  with  that  of  an  infinite  array,  except  near 
endfire.  For  this  reason,  we  often  use  the  65  x  149  array  results  instead  of  the  corresponding 
(and  nearly  Identical)  infinite  array  results. 

Nonisolating  feeds  can  be  represented  by  an  (N  +  M)  x  (N  +  M)  scattering  matrix  [Sp],  where 
N  is  the  number  of  antenna  ports  and  M  is  the  number  of  input  ports.  This  scattering  matrix  can 
be  written  in  partitioned  form  as 


isp] 


F 

o  F 

11 

j  S12 
| _ 

F 

!  sf 

21 

\  “22 

(27) 


where  §  .  is  an  M  x  M  matrix  which  represents  the  input  mismatches  and  cross 'Couplings  among 
the  input  lines,  §21  is  an  N  x  M  matrix  representing  the  coupling  from  the  input  ports  to  the  an¬ 
tenna  ports  (the  antenna  illumination  under  matched  conditions),  and  is  an  N  x  N  matrix  giving 

the  output  mismatches  (looking  back  at  the  feed)  and  the  output  cross -couplings  among  the  feed 

F 

lines.  When  only  reciprocal  devices  are  used  in  constructing  the  feed,  S,2  will  be  just  the  trans- 
pose  of  fi21 .  Also,  the  feed  scattering  matrix  will  be  unitary  (the  inverse  of  [Sp)  is  given  by  the 
transpose-conjugate  of  [Sp]),  if  the  components  in  the  feed  are  lossless. 

Thus,  for  a  particular  feed  network,  the  characteristics  of  the  components  which  comprise 
the  network  can  be  combined  to  obtain  a  scattering  matrix  in  the  form  of  Eq.  (27).  Specific  ex¬ 
amples  are  treated  in  Sec.  III-C.  in  general,  we  can  state  that  the  performance  of  these  non- 
isolating  feed  systems  is  quite  sensitive  to  the  design  parameters  of  the  network  (principally  line 
lengths)  because  of  multiple  reflections  within  the  structure.  These  reflection-  can  add  vectori- 
ally  at  the  various  ports  to  give  either  very  good  or  very  poor  array  illumination  depending  on  the 
particular  design  of  the  network. 
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Fig.  23.  Typical  array  function  ihowing  spurious  lobss  for  |  |  ■  0. 1 

for  an  array  ut'ng  throe-port  junctions. 

The  phase  shifters  used  in  the  composite  system  (antennas,  feed  network,  and  phase  shifters) 
can  also  have  a  strong  influence  on  the  array  illumination.  For  example,  when  reciprocal  phase 
shifters  are  used,  the  part  of  the  wave  reflected  from  an  antenna  that  is  reflected  back  to  the  an¬ 
tenna  ty  the  feed  network  generates  a  spurious  beam  in  a  different  direction,  resulting  in  a  high 
sidelobe  in  the  array  pattern,16  as  indicated  in  Fig.  23.  Quantitatively,  for  any  feed  with  reflec¬ 
tion  coefficient  rA  looking  into  the  antennas  and  looking  into  the  feed  (see  Fig.  21),  the  illu¬ 
mination  coefficients  are  related  to  the  original  coefficients  an  exp[-jnkD  sine)  by 
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(28) 


an  =  an  2  <-4)q  (rArB>q  exp(-jnkEK2q  +  1)  sinG] 

q  =  0 

Thus,  for  a  reflection  coefficient  product  of  magnitude  0.1,  the  resultant  far-field  pattern  would 
be  approximately  as  shown  in  Fig.  23  for  a  phase  shift  setting  corresponding  to  a  three -beam width 
pointing  angle  from  broadside.  On  the  other  hand,  when  nonreciprocal  phase  shifters  are  used, 
the  re -reflected  waves  arrive  back  at  the  antennas  with  the  same  relative  phase  shifts  as  the  pri¬ 
mary  waves  and  hence,  to  a  first  order,  only  affect  the  amplitude  (gain)  of  the  array  pattern. 


B.  Arrays  Fed  by  Independent  Sources 

In  this  section,  we  examine  some  small-array  configurations  to  determine  (1)  the  correlation 
of  finite  and  infinite  array  element  impedance  and  element  gain  properties  and  (2)  the  magnitude 
and  extent  of  edge  effects  in  small  arrays. 

1.  Impedance  Behavior  of  Linearly  Polarized  Arrays 


For  these  investigations,  we  will  restrict  ourselves  to  7  x  9  and  9X11  arrays.  (The  im¬ 
pedance  matrix  for  a  9  x  1 1  array  is  the  largest  which  can  be  conveniently  handled  within  the  core 
storage  of  an  IBM  7094  digital  computer.)  The  geometry  and  element  numbering  system  we  will 
adopt  for  a  9x11  rectangular -grid  array  are  shown  in  Fig.  24.  Alternate  rows  of  the  array  are 
offset  by  D */Z  to  obtain  the  triangular-grid  array  discussed  later. 

The  impedance  variation  with  scan  angle  for  the  center  element  of  a  7  x  9  array  of  x/2  dipole 

radiators  above  a  ground  plane  (D  =  D  =  0.5X,  s  =  0.25A)  is  shown  in  Fig.  25.  One  set  of  curves 

x  y 


Fig.  24.  Geometry  and  element  numbering. syitem  Fig.  25.  Comparison  of  impedance  variation  of  the 

for  a  9  X  11  element  array  configuration  (element  center  element  In  a  7  X  9  array  of  V?  dipoles  with 

50  Is  the  center  element  for  this  array).  that  of  on  infinite  array. 
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gives  the  results  for  the  exact,  matrix  inversion  solution,  wliile  another  set  gives  the  results 
obtained  with  the  approximation  of  Eq.  (26).  Both  sets  of  curves  are  then  compared  with  the 
Impedance  variation  of  a  large  array  (6$  x  149)  with  the  same  geometry.  It  is  evident  that  the 
7X9  array  results  (both  exact  and  approximate)  agree  reasonably  well  with  the  large -array 
results  in  all  three  scan  planes  (the  difference  in  the  magnitudes  of  the  7X9  array  impedance 
and  the  large-array  impedance  at  a  particular  scan  angle  is  less  than  14  percent).  Near  broad¬ 
side,  the  small-array  impedances  exhibit  rather  unusual  behavior,  while  at  larger  scan  angles, 
they  tend  to  oscillate  around  the  values  for  a  large  array.  An  increase  In  array  size  to  9  x  11 
results  in  a  somewhat  better  approximation  to  an  infinite  array,  but  the  deviations  in  the  magni¬ 
tude  of  the  impedance  compared  to  the  large  array  are  still  of  the  order  of  10  percent.  A  still 

better  approximation  to  the  infinite  array  impedance  variation  can  be  found  by  taking  the  geo- 

17 

metric  mean  of  the  7X9  and  9  x  11  impedances  at  each  of  the  scan  angles. 

Similar  results  are  observed  when  the  arrays  have  triangular  grids  as  shown  in  Fig.  26.  In 
this  figure,  the  impedance  variation  with  scan  angle  of  a  central  element  of  a  9  x  11  equilateral 
triangular-grid  array  (D  =  0.6  X,  s  =  0.25  X)  is  compared  to  that  of  a  corresponding  infinite  array 
for  three  scan  planes.  Again  it  is  found  that  the  error  in  approximating  the  impedance  behavior 
of  a  large  array  by  a  small  array  is  only  about  10  percent.  The  small -array  impedance  varia¬ 
tion  oscillates  around  that  of  an  infinite  array  which  is  again  just  typical  of  the  behavior  expected 
when  a  function  is  approximated  by  a  truncated  Fourier  series. 

The  small-array  behavior  differs  more  markedly  from  that  of  an  infinite  array  when  the 
ground  plane  is  removed,  as  evidenced  by  Fig.  27  where  the  impedance  variation  of  the  center 
element  in  7  x  9  and  9X11  arrays  is  compared.  It  is  clear  from  the  very  poor  correlation 


Fig.  26.  Comparison  of  impedance  variation  of  the 
center  element  in  a  9X11  equilateral  triangular- 
grid  array  of  >/2  dipole*  with  that  of  an  Infinite 
array. 


Fig.  27.  Comparison  of  impedance  variation  of  the 
center  element  in  7  X  9  and  9x11  square -grid  ar¬ 
ray*  of  V?  dipoles  without  ground  plane*. 
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between  the  two  array  impedance  variations  that  small  arrays  of  dipoles  without  a  ground  plane 
or  slots  in  a  ground  plane  are  not  very  useful  for  predicting  the  performance  of  large  arrays. 

As  discussed  above,  the  impedance  variation  (magnitude)  with  scan  angle  for  the  center  ele¬ 
ment  (No.  50  of  Fig.  24)  of  a  9  x  H  square-grid  array  agrees  within  about  10  percent  with  the 
values  for  an  infinite  array.  For  the  impedance  variation  with  scan  angle  of  an  interior,  non¬ 
center  element,  there  is  still  reasonably  good  agreement  with  an  infinite  array,  even  though 
there  may  be  substantial  asymmetry  in  a  given  scan  plane:  for  a  particular  scan  plane  speci¬ 
fied  by  <p  =  0 the  impedance  variation  with  scan  angle  0  for  one  direction  of  scan  from  broad¬ 
side  will  differ  from  that  when  the  beam  is  scanned  in  the  opposite  direction  ZD  (<p.  +  *■,  6)  ^ 

ZD  «VG>- 
m 

The  situation  is  considerably  worse  when  the  elements  are  on  the  edges  of  the  array  as 
shown  in  Fig.  28  for  the  element  (46)  in  the  middle  of  the  edge  where  the  dipoles  are  parallel,  and 
in  Fig.  29  for  the  element  (5)  in  the  center  of  the  edge  where  the  dipoles  are  collinear.  In  these 
figures,  it  is  observed  that  the  correlation  of  the  impedance  variation  with  that  for  an  infinite  ar¬ 
ray  is  very  poor.  In  addition,  there  is  considerable  asymmetry  as  evidenced  by  the  fact  that 

ZD  (<f>.  +  ir,  0)  differs  markedly  (in  magnitude  and  phase)  from  ZD  ( <f> .,  0).  Note  that  the  array 
um  um 

is  symmetrical  about  element  46  for  H-plane  scan,  and  therefore  no  asymmetry  is  observed  in 

the  H-plane  impedance  variation  with  scan  angle,  as  shown  in  Fig.  28.  Similarly,  the  array  is 

symmetrical  about  element  5  for  E-plane  scan. 

From  these  results,  we  conclude  that  for  a  9  x  11  (or  smaller)  array,  there  is,  strictly 
speaking,  no  typical  element;  that  is,  every  element  has  a  different  impedance  variation  with 
scan  angle.  However,  for  dipoles  above  ground,  the  center  element  and  the  edge  effects  com¬ 
puted  or  measured  for  the  small  array  would  give  a  reasonable  indication  of  those  to  be  expected 
in  a  large  array.  For  dipoles  without  a  ground  plane,  slots  in  a  ground  plane  or  open  waveguides, 
the  9xll  array  gives  only  rough  qualitative  estimates  of  large-array  behavior. 


2.  Pattern  Behavior  of  Regular  Linearly  Polarized  Arrays 

Since  we  can  model  the  excitation  of  a  typical  element  of  an  independently  fed  array  by  a 
circuit  such  as  that  of  Fig.  22,  we  can  conceptually  turn  on  each  generator  of  the  array  individ¬ 
ually,  with  the  others  turned  off  (e^  =  0,  m  ^  n).  When  the  generators  are  shut  off,  all  passive 
elements  are  effectively  terminated  in  an  impedance  Z^.  With  only  eR  ±  0,  a  current  in  will 
flow  past  some  reference  point  on  antenna  n,  directly  exciting  antenna  n  and  parasitically  ex¬ 
citing  the  rest  of  the  array  through  mutual  coupling.  A  measurement  of  the  far  field  under  these 

conditions  produces  some  relative  field  strength  which  we  will  call  f  ($,  9),  since  it  is  due  to 

th  n 

the  n  element  generator.  The  experiment  can  be  repeated  for  each  element,  and  the  far  field 

can  be  found  by  taking  the  vector  sum  of  the  individual  element  contributions,  properly  weighted 

by  a  phase  term  dependent  on  the  element  location 


F(*.  0)  »  £  in  e*PlJk?n  *  *R)  (29) 

n 

where  SD  is  a  unit  vector  in  the  direction  of  the  observation  point,  and  5  is  a  vector  from  the 
origin  to  the  n  element  location.  For  completely  arbitrary  distributions  of  elements,  little 
more  can  be  said.  However,  for  the  important  case  in  which  (a)  the  array  is  large  compared  to 
the  region  over  which  element  interaction  is  strong  so  that  edge  effects  can  be  neglected,  (b)  the 
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Fig.  28.  Comparison  of  imp* dance  variation  of  the 
element  in  the  center  of  the  parallel  dipole  edge  of 
a  9  X  l!  array  of  \/2  dipoles  with  that  of  a  large 
array. 


Fig.  29.  Comparison  of  impedance  variation  of  the  element  in  the  center  of  the  col  I  inear 
dipole  edge  of  o  9  X  1 1  array  of  dipoles  with  that  of  o  large  array. 


elements  and  their  generator  impedances  are  nominally  identical,  and  (c)  the  elements  are  reg¬ 
ularly  spaced  on  a  flat^  surface  so  that  all  interior  elements  "see"  the  same  interaction  environ¬ 
ment,  the  pattern  of  each  element,  measured  ao  indicated  above,  will  be  essentially  the  same  as 
the  pattern  of  any  other  element.  We  can  then  factor  Eq.  (29)  and  write 

F (<#>,  0)  =  f(4>,  G)  £  in  exp  [jkpn  •  eR]  (30) 

n 

where  f (p,  0)  is  a  typical  f ^(0,  9). 

The  factor  f{0,  0)  multiplying  the  summation  in  Eq.  (30)  is  the  pattern  of  a  typical  (central) 
element  in  the  presence  of  all  the  other  elements  when  they  are  terminated  in  the  impedance  from 
which  they  are  nominally  excited,  and  is  usually  referred  to  as  the  "element  factor."  The  i  of 
Eq.  (30)  represents  the  current  flow  in  the  n  n  element  due  solely  to  the  excitation  of  its  own  gen¬ 
erator.  For  all  elements  (neglecting  edge  effects),  i  is  related  to  en  of  Fig.  22  by  the  same  pro¬ 
portionality  independent  of  the  other  element  excitations,  and  we  could  replace  i  by  en  in  Eq.  (30), 
except  for  a  scale  factor  which  is  the  passive  element  driving  impedance. 

The  practical  significance  of  this  development  lies  in  the  fact  that  the  summation  of  Eq.  (30) 
is  exactly  the  array  factor  of  an  array  of  isotropic  radiators,  the  properties  of  which  have  been 
extensively  discussed.  The  element  factor  contains  essentially  all  the  pattern  effects  of  the  ele¬ 
ment  type  ind  the  interaction  between  elements.  The  pattern  of  the  entire  array  when  only  one 
element  is  excited  and  the  others  terminated  is  much  broader  than  the  pattern  when  the  entire  ar¬ 
ray  is  excited;  hence,  f(<p,  0)  is  a  function  whose  variation  with  angle  is  slow  compared  to  any 
reasonably  directive  array  factor  in  large  arrays.  The  patterns  are  qualitatively  similar  to  those 
of  Fig.  30.  The  element  factor  acts  as  a  window  of  varying  degrees  of  opaqueness  through  which 
the  array  factor  "looks."  We  can  qualitatively  conclude  that  under  normal  circumstances,  in  large 
arrays,  the  element  iVctor  structure  will  have  little  effect  on  the  relative  structure  of  the  main 
beam  and  close-in  sidelobes  of  a  large  array.  However,  in  arrays  in  which  the  beam  is  scanned 
by  element  phasing,  the  absolute  strength  of  the  pattern  will  be  varied  with  angle  to  conform  to 
the  element  factor  weighting,  since  the  scanning  moves  the  array  factor  in  space  while  the  element 
factor  is  stationary,  as  indicated  by  Fig.  30. 

For  elements  in  which  mutual  impedances  can  be  analytically  determined,  the  computation 
of  the  element  factor  is  straightforward.  Perhaps  more  important,  in  the  frequent  case  of  ele¬ 
ments  for  which  computational  formulas  do  not  exist,  the  element  factor  can  be  determined  ex¬ 
perimentally  by  building  an  array  only  large  enough  to  render  edge  effects  negligible  on  the  pat¬ 
tern  of  the  central  element.  Finally,  we  can  also  rationalize  some  generally  valid,  quantitative 
conclusions  about  element  factor  shapes. 

Of  course,  if  the  antenna  elements  are  separated  far  enough  so  that  the  interactions  are  neg¬ 
ligible,  the  fn($,0)  become  the  patterns  of  isolated  elements  (which  are  given  in  the  literature 
for  a  wide  variety  of  radiators).  However,  only  when  the  elements  are  separated  by  distances 
that  are  much  greater  than  their  largest  dimension  is  it  safe  for  us  to  assume  that  interaction 
effects  are  negligible.  Unfortunately,  achieving  freedom  from  interaction  in  this  way  results  in 
spacing  the  elements  sufficiently  far  apart  that  several  grating  lobes  of  the  array  factor  appear 
in  visible  space,  as  indicated  in  Sec.II-D. 

t  The  modification  for  large,  regularly  (paced  array*  on  curved  (urfoce*  with  large  radii  of  curvature  (compared 
to  X)  i*  tedious,  but  not  conceptually  difficult 
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Fig.  30.  Qualitative  relationship  between  element  pain  function  and  array  factor  ({  =  angle  from 
broadside):  (a)  array  factor,  (b)  element  gain  function,  (c)  resulting  pattern. 

The  more  interesting  case  for  many  array  applications  is  that  of  closely  spaced  elements  - 
in  particular,  the  case  of  only  slightly  directive  elements  spaced  less  than  a  wavelength  on  cen¬ 
ters.  Here  some  element  factor  shapes  can  be  inferred  from  a  consideration  of  the  gain  and 
directivity  variation  of  electronically  scanned  arrays. 

The  array  gain  G  is  a  measure  of  efficiency  of  the  array  in  using  all  the  energy  available, 
and  involves  questions  of  array  impedance  matching  and  ohmic  losses.  Although  mismatch  ef¬ 
fects  are  not  always  included  as  a  factor  in  the  gain  of  an  antenna,  they  will  be  included  in  this 
discussion.  Such  effects  are  taken  into  account  by  defining  the  gain  as  the  ratio  of  power  density 
per  unit  solid  angle  at  the  peak  of  the  main  beam  to  the  power  density  that  would  be  achieved  by 
radiating  all  the  available  power  isotropically.^  For  a  large,  flat,  equally  spaced  array  of  iden¬ 
tical  radiators,  it  is  a  simple  matter  to  determine  the  gain  as  a  function  of  scan  angle  in  terms 
of  the  element  factor  shape. 

The  power  density  per  unit  solid  angle  is 

PIVV  |B|2|,l*o’eo»l2  Zen 

n 

where  B  is  an  angle  invariant  factor,  the  value  of  which  will  be  seen  to  be  immaterial.  The 
power  available  from  each  generator,  assuming  that  the  equivalent  circuit  of  Fig.  22  applies,  is 

tThe  IEEE  Standard*  define*  gain  a*  relative  to  the  total  power  into  the  antenna.  Ignoring  unmatched  effect*. 

The  definition  of  gain  used  herein  thu*  differ*  from  the  itondord  definition  by  the  factor  I  -  |  T  |  2, 
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I  e  I  /4R  ,  where  R  =  R  (Z  ),  and  we  assume  all  element  drives  have  the  same  impedance. 

1  n'  '  g  g  e  g  2 

Hence,  the  total  power  available  to  the  array  is  £  je  |  /4R  from  which  it  follows  that  the  gain 

n  n  ® 

of  an  array  is 

l£  en|2 

O(*0.eoi  =  (4»)  |B|z4RgUl»0,e0ll  p  ■ 

Note  that  if  we  drive  only  one  element  in  the  array,  the  gain  of  that  element  is  given  by 
g(<P0.00)  =  (4*)4Rg|B|2  |f(0o,0o)|2  . 

Further,  if  =  i,  for  all  n,  then 


*  e. 


n 


n' 


2knl' 

n 


=  N 


where  N  is  the  number  of  elements  in  the  array.  For  arbitrary  taper,  we  can  define  a  taper 
efficiency  tj  by 


£ 

n 


2 

7 


TJ  N 


Comparison  of  the  relationships  for  G,  g,  and  rj  N  gives  the  simple  result 

G(«,9)  =  g(«f>,Q)  rjN  .  (31) 

The  factor  g(tp,  0)  is  referred  to  as  the  "clement  gain  function,"  and  its  spatial  variation  with  the 
beam -pointing  angles  <f>  and  €*  is  |f($,0)|2.  Equation  (31)  is  the  quantitative  statement  of  the 
fact  that  the  element  factor  acts  as  a  weighting  function  on  the  field  strength  of  the  array.  It  is 
easily  established  that  the  directi\ity  of  a  planar  array  of  N&  active  radiators  with  no  visible 
grating  lobes  is 

Na 

U(0)  =  4irrj  — =•  b  cos  0  (32) 

A4 


where  0  is  the  angle  from  the  array  broadside  direction,  and  a  is  the  area  allotted  to  each  ele¬ 
ment.  For  independently  excited  ri  diators  with  no  ohmic  loss,  which  can  support  only  a  single 
polarization,  the  array  directivity  and  gain  as  defined  herein  can  differ  only  as  the  result  of  mis¬ 
match  losses.  These  losses  can  be  taken  into  account,  in  terms  of  the  voltage  reflection  coeffi¬ 
cient  r(<*>,0)  seen  when  looking  into  a  typical  element  when  the  entire  array  is  excited,  by  writing 


g<*, ©)  =  u<e)  (i  -  |r(*,e),?'  . 

Substitution  of  Eq.  (32)  into  this  expression  sheds  some  light  on  the  element  factor  behavior  to 
be  expected  as  a  function  of  element  spacing: 


g(*,e)  =  4*  A  cose  (i-|r(*,e)|2) 

A* 


(33) 
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Fi0.  31 .  H-plane  element  gain  functions  for  ari 


We  can  now  draw  some  inferences  about  the  general  nature  of  the  gain  function  shapes  to  be  ex- 

t 

pected  in  planar  arrays. 

The  only  place  the  type  of  radiator  enters  the  above  expression  is  implicitly  in  r(<£,  8).  Re¬ 
gardless  of  the  type  of  radiator,  we  can  choose  the  generator  impedance  to  match  the  element 
at  some  angle  of  scan  of  our  choice.  If,  for  the  sake  of  illustration  (and  maximum  broadside 
gain),  we  match  the  elements  when  the  array  is  phased  for  a  broadside  beam,  we  can  make  the 
following  generalization  about  the  element  pattern  for  small  (less  than  a  wavelength  in  dimension) 
elements  from  Eq.  (33).  As  we  space  the  elements  farther  apart,  the  value  of  g(0,  0)  can  be  made 
to  increase  directly  with  the  area  per  elemeni  if  the  proper  generator  impedance  is  chosen  for 
each  spacing.  However,  as  the  element  spacing  is  increased,  the  angles  at  which  grating  lobes 
become  visible  decrease.  At  these  angles,  the  array  directivity  must  drop  rapidly  by  approxi¬ 
mately  3  db,  and  the  gain  function  must  behave  in  a  similar  fashion  (we  are  not  matched  at  these 
angles,  so  only  approximate  statements  are  justified).  That  the  element  pattern  does  behave 
in  such  a  manner  has  been  verified  analytically  for  some  elements  and  experimentally  for  many 
others.  Figure  31  shows  element  patterns  of  the  center  element  of  an  array  of  dipoles  without 
a  ground  plane  for  various  spacings  on  a  square  element  grid,  and  Fig.  32  shows  similar  data  for 
dipoles  0.25  X  above  a  ground  plane.  The  figures  illustrate  the  effects  stated. 

If  the  elements  are  not  matched  to  maximize  g(0,  0)  for  each  spacing,  the  result  is  a  decrease 
in  broadside  gain  and,  usually,  a  slight  increase  in  the  width  of  the  element  factor  main  lobe.  For 
example,  a  convenient  generator  impedance  choice  from  a  practical  standpoint  is  to  match  the  im¬ 
pedance  of  a  single  isolated  antenna  (with  ground  plane,  if  used).  Indeed,  Eq.  (5)  indicates  that 
for  an  essentially  infinite  array  scanned  over  all  phasings  (including  those  producing  grating 
lobes),  the  average  driving  impedance  is  just 

ir 

j'j'zD(«.fl)d^  =  Z00  00  134) 

-IT 

which  is  essentially  the  same  as  the  impedance  of  an  isolated  element  for  many  practical  types. 
Consequently,  this  choice  comes  close  to  minimizing  the  maximum  impedance  variation,  and  will 
usually  produce  a  match  at  some  scan  angle  other  than  broadside.  The  resulting  element  gain 
function  is  usually  "saddle-shaped"  as  indicated  in  Fig.  33, 

For  arrays  on  curved  surfaces,  if  the  curvature  is  shallow  compared  to  the  extent  over  which 
appreciable  coupling  takes  place,  the  element  gain  function  shape  relative  to  each  element's  own 
broadside  angle  will  be  similar  to  those  of  an  equivalently  spaced  planar  array.  Thus,  even 
though  Eq.  (30)  fails  to  apply,  the  element  patterns  can  usually  be  taken  into  account  in  a  straight¬ 
forward  and  relatively  simple  manner. 

The  element  gain  function  concept  has  both  practical  and  fundamental  significance.  It  is  a 
practical  way  to  assess  the  gain-vs -scan -angle  performance  of  a  large  array  by  using  a  small 
or  modest  array  without  the  requirement  of  a  complex  feed  system.  We  need  only  construct  an 
array  large  enough  so  that  the  central  element  "sees"  negligible  edge  effects  (how  large  in  each 
dimension  depends  on  the  rate  of  decay  of  coupling  for  the  particular  element  in  the  direction  in 


t  Equation  (31)  applies  also  to  linear  arrays;  hence,  a  dependence  analogous  to  Eq.  (33)  of  the  gain  function  of 
an  element  in  a  linear  array  on  element  spacing  can  be  made  (see  Fig. 43),  but  the  result  is  complicated  by  the 
fact  that  the  variation  of  array  gain  with  scan  angle  depends  on  the  element  pattern  in  the  plane  orthogonal  to 
the  array. 
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Fig.  33.  Principal  plane  gain  function*  for  a  'S/I  dipole  in 
an  array;  generator  matched  to  the  impedance  of  an  isolated 
element. 


Fig.  34.  Hybrid-fed,  cro**ed-dipo!e-pair  element  configuration;  typically,  Z^  =  -  Z^. 
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question,  as  elaborated  on  previously).  All  elements  are  then  terminated  except  the  central 
element,  and  a  pattern  i3  taken.  Experiments  in  design  parameters  can  be  conducted  on  a  lim¬ 
ited  model. 

There  is  also  a  connection  between  the  shape  of  the  element  gain  function  and  the  complex 

18 

driving  impedance  of  the  elements.  For  example,  note  that  in  an  essentially  infinite  array  for 
which  the  elements  have  a  voltage  retiection  coefficient  behavior  T (0,0),  the  current  into  an 
element  at  scan  angle  <t>,  0  is  1  -  T (<p,  0).  Furthermore,  the  shape  of  the  element  gain  function 
can  be  used  to  demonstrate  a  fundamental  limitation  on  the  average  reflected  power  over  all 
phasings,  applicable  to  any  element.1 

3.  Gain,  Polarization,  and  Impedance  Properties  of  Independently  Fed 
Circularly  Polarized  Arrays 

An  additional  effect  of  element  interaction  on  array  patterns  is  found  in  arrays  of  elements 
which  can  support  more  than  one  polarization.  The  general  effect  of  the  coupling  is  to  cause  the 
array  polarization  to  vary  with  scan  angle  in  a  manner  dependent  on  all  the  factors  that  enter  in¬ 
to  the  determination  of  coupling  effects:  element  type,  feed  circuit  details,  and  element  grid 
shape  and  spacing. 

To  illustrate  these  effects,  we  present  in  this  section  some  results  of  an  analysis  of  arrays 
of  elements  made  of  crossed  dipoles.  In  particular,  we  examine  arrays  of  such  elements  driven 
through  a  hybrid  power  divider  as  indicated  in  Fig.  34,  so  that  the  elements  can  be  considered 
independently  excited.  In  such  a  case,  an  element  gain  function  for  a  dipole  pair  in  the  array 
can  be  defined  in  the  same  manner  as  for  linearly  polarized  arrays  in  the  previous  section.  How¬ 
ever,  there  are  contributions  to  the  gain  function  from  both  co-directed  and  orthogonally  directed 
dipoles,  and  hence  the  gain  function  must  be  defined  relative  to  some  polarization  vector  ep 

g(*.  ©,  e^)  =  JT  I  (ip  ‘  fi‘(*'  0)  {1  ~  r*(*'  0)5 

-j(e^  •  ep  fy(*.  0)  (1  -  Ty(0,  0))  |  2  (35) 

for  an  interior  element  of  a  large  array,  where  f.xe“  is  the  vector  element  pattern  of  the  x- 

x  *  x 

directed  elements,  and  T  (0,0)  is  the  reflection  coefficient  looking  into  those  elements,  and 
similarly  for  the  y-directed  elements.  It  is  readily  shown  that  for  an  infinite,  regularly  spaced 
array  above  a  ground  plane,  the  broadside  gain  relative  to  a  right-circular  sampling  polariza¬ 
tion  vector  is 

g(0.  0,  e”*  )  =  4ir  -% 
r  A 

when  all  dipoles  are  matched  to  the  feed  hybrid  at  broadside  (Ry  =  R^JO,  0)  =  RjpO,  0)),  a  result 
expected  on  the  basis  of  directivity  considerations. 

From  Eq.  (35),  it  is  possible  to  derive  expressions  for  the  gain  relative  to  right -circular, 
left-circular,  or  any  arbitrary  linear  sampling  polarization  vector  at  any  scan  angle.  These 
relations  can  then  be  used  to  determine  the  characteristics  (ellipticity,  tilt  angle,  left-  and 
right  -circular  gain,  etc.)  of  the  wave  transmitted  in  any  direction. 

In  order  to  assess  the  effects  of  an  array  environment  on  the  dipole  gain  and  impedance,  we 
first  consider  an  isolated  dipole  pair.  The  polarization  of  the  transmitted  wave  at  various  points 
in  space  is  the  characteristic  of  importance  for  the  isolated  element.  A  typical  power  pattern  as 
would  be  measured  by  a  spinning,  linearly  polarized  receiving  dipole  is  shown  in  Fig.  35  for  a 
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dipole  pair  above  an  infinite  ground  plane.  Also  shown  are  the  maximum  and  minimum  bounds 
of  the  corresponding  pattern  for  a  crossed-dipole  pair  at  any  point  in  space. 

When  seven  elements  are  arranged  in  a  linear  array  (no  ground  plane),  the  polarization 
characteristic  of  the  gain  function  for  the  center  element  is  as  shown  in  Fig.  36  for  a  cut  in  the 
direction  of  the  array  when  the  dipoles  are  matched  to  the  feed  circuit  at  broadside  and  the  fourth 
port  of  the  hybrid  (Fig.  34)  is  terminated  in  a  matched  load.  Note  that  the  array  environment  has 
substantially  changed  the  gain  from  that  of  an  isolated  element,  with  the  result  that  nearly  circu¬ 
lar  polarization  is  obtained  at  a  scan  angle  of  45°,  while  at  broadside  the  polarization  ellipse 

axial  ratio  (max.  gain/min.  gain)  is  about  3db.  (This  is  similar  to  the  results  observed  by 
1 9 

Parad  and  Kreutel  for  a  seven -element  hexagonal  array  of  circular  waveguides.)  The  polari¬ 
zation  characteristic  can  be  improved  only  slightly  by  rotating  the  elements  45°  as  demonstrated 
by  the  dashed  curves  in  Fig.  36.  Note  that  the  polarization  of  the  transmitted  wave  at  broadside 
is  elliptical  even  with  the  dipoles  rotated. 

As  a  final  result,  we  consider  the  center  element  of  a  13  x  13  square-grid  array  above  a 
ground  plane.  A  principal  plane  (<p  =  0°)  gain  pattern  as  would  be  measured  by  a  spinning,  re¬ 
ceiving  dipole  when  the  elements  are  matched  for  broadside  radiation  and  the  fourth  port  is  ter¬ 
minated  is  shown  in  Fig.  37.  For  this  array,  perfect,  right -circular  polarization  is  obtained  at 
broadside  even  in  the  presence  of  mutual  coupling.  This  is  a  consequence  of  the  complete  sym¬ 
metry  of  the  array  about  the  center  element.  The  polarization  properties  in  other  planes  of 
scan  for  this  array  are  shown  in  Fig.  38,  where  the  polarization  ellipse  axial  ratio  variation  with 
scan  angle  for  both  an  isolated  dipole  pair  and  the  center  dipole  pair  in  the  13  x  13  array  is  plotted. 
Note  that  the  array  environment  can  substantially  alter  the  polarization  of  the  transmitted  wave 
(as  compared  to  that  of  an  isolated  dipole  pair)  and  can  in  fact  improve  the  circularity  of  the  ra¬ 
diated  wave  for  optimum  ground  plane  spacing  (s  »  0.25  X  for  D  =  D  =  0.6  X).  Since  the  center 

x  y 

element  of  the  13  x  13  array  "sees"  an  environment  which  is  approximately  the  same  as  that  seen 
by  an  element  in  an  infinite  array,  the  polarization  results  given  in  Fig.  38  should  be  quite  indic¬ 
ative  of  large-array  polarization  variation  with  scan  angle  to  be  expected  from  an  array  of 


Fig.  38.  Polarization  allipM  axial  ratio*  (a)  holatad  clamant,  (b)  cantor  clamant  in  13 X  13  array. 
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Fig.  39.  Input  port  impodanco  variation  for  the  cantor  a! ament  of  a  13  X  13  array 
of  hybrid-fad,  crowed -dipole  pain. 
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cross -polarized  elements  fed  by  a  hybrid.  In  addition,  the  input  impedance  variation  with  scan 
angle  at  port  i  of  Fig.  34  and  the  power  dissipated  (relative  to  the  power  available  from  the 
source)  in  the  orthogonal  port  termination,  as  shown  in  Figs.  39  and  40,  should  be  within  approx¬ 
imately  10  percent  of  that  which  would  be  observed  in  a  very  large  array. 

From  the  above  results,  it  is  evident  that  mutual  coupling  effects  can  substantially  alter  the 
characteristics  (gain,  impedance,  and  polarization  variation  with  scan  angle)  of  a  circularly  po¬ 
larized  element  when  it  is  placed  in  an  array  environment.  However,  it  is  also  clear  that  a  prop¬ 
erly  designed  array  of  hybrid-fed,  crossed-dipole  pairs  will  give  good  performance  for  moder¬ 
ate  scan  angles  (0  *  40*).  Observe  that  the  characteristics  of  an  element  in  an  array  depend 

max 

on  the  array  geometry  for  both  linearly  polarized  and  circularly  polarized  dipole  arrays;  for  ex¬ 
ample,  there  is  an  optimum  height  above  the  ground  plane  for  every  inter-element  spacing.  It 
is  also  worth  noting  that  the  input  impedance  as  seen  by  a  source  driving  an  element  in  an  array 
varies  somewhat  less  with  scan  angle  for  crossed-dipole  arrays  (compared  to  linearly  polarized 
arrays),  particularly  in  the  reactive  part. 

It  should  be  stressed  that  these  results  are  dependent  on  the  circuit  used  toobtain  the  two  or¬ 
thogonal  excitations.  For  example,  if  the  fourth  port  of  the  hybrid  were  left  open-circuited,  the 
power  dissipated  in  that  port  when  it  is  terminated  (Fig.  40)  would  be  reradiated  instead  of  absorbed. 
Examination  of  the  phases  involved  shows  that  such  radiation  is  of  the  opposite  sense  circular  po¬ 
larization  to  that  desired,  resulting  in  greater  depolarization  than  indicated  by  Figs.  36-38.  In 

jq 

fact,  it  has  been  pointed  out  '  that  to  minimize  depolarization  effects  due  to  coupling  in  more  than 
one  polarization,  it  is  necessary  to  provide  a  matched  load  to  all  polarizations,  even  if  only  one 
is  desired.  For  example,  in  an  array  designed  solely  to  radiate  one  sense  of  circular  polari¬ 
zation,  a  load  must  be  provided  to  the  other  sense  to  minimize  depolarization. 

4.  Pattern  Behavior  of  Irregular  Arrays 

The  possibility  of  shaping  the  far-field  patterns  of  arrays  with  equal  amplitude  excitation 
by  unequally  spacing  the  elements  has  been  extensively  studied,  but  almost  exclusively  under  the 
assumption  that  the  effect  of  the  mutual  coupling  between  the  elements  can  be  ignored.^  There 
is  experimental  and  theoretical  evidence,2®  however,  that  the  effect  of  coupling  on  the  patterns 
of  such  arrays  is  not  generally  negligible.  For  example.  Fig.  41  compares  the  pattern  of  an 
unequally  spaced  array  of  Galejs21  as  computed  ignoring  coupling,  and  as  computed  assuming 
the  array  to  be  constructed  of  parallel  short  dipoles  0.25  \  above  ground,  independently  fed  by 
sources  of  the  type  shown  in  Fig.  22,  with  chosen  to  match  the  impedance  of  an  isolated  dipole 
0.25  A  above  ground.  The  patterns  shown  are  for  a  broadside  beam;  the  relative  pattern  (most 
notably,  the  near-in  sidelobe  ratio)  also  varies  with  scan  angle  in  the  presence  of  coupling  as 
indicated  by  the  corresponding  patterns  of  Fig.  42. 

The  necessity  for  including  coupling  in  pattern  computations  for  unequally  spaced  arrays  is 
in  marked  contrast  to  the  fact  that  the  effects  of  coupling  can  safely  be  ignored  in  relative  pat¬ 
tern  computations  for  large  equally  spaced  arrays.  This  latter  simplification  is,  of  course,  due 
to  the  fact  that  for  regular  arrays  large  enough  to  ignore  edge  effects,  all  elements  in  the  array 
"see"  the  same  coupling  environment  and,  hence,  have  the  same  input  impedance. 

TK*re  are  only  two  cases  in  which  it  is  certain  that  the  effects  of  such  coupling  on  the  shape 
of  the  pattern  are  negligible.  First,  these  effects  can  certainly  be  ignored  when  the  elements 

t  See  Ref.  2)  for  a  recent  bibliography. 
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Fig.  41.  Comparison  of  computed  parallel  dipole  array  factors  for  an  unequally  spaced 
array;  broadside  beam. 


Fig.  42.  Comparison  of  computed  parallel  dipole  array  factors  for  an  unequally  spaced 
array;  45*  scan  angle  frcm  broadside  beam. 
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are  fed  by  ideal  constant -current  generators  so  that  their  excitation  is  independent  of  their  im¬ 
pedance.  However,  such  generators  are  hard  to  find.  Second,  coupling  can  obviously  be  ignored 
when  essentially  all  elements  of  an  array  are  widely  separated  compared  to  the  aperture  area  of 
such  r.n  element.  This  is  not  usually  an  acceptable  solution  in  practical  arrays,  since  unequal 
spacing  allows  control  of  the  sidelobes  commensurate  with  that  of  an  equally  spaced,  amplitude 
tapered  array  only  over  a  limited  angle  about  the  main  beam.  This  region  is  out  to  an  angle 
approximately  equal  to  that  at  which  a  grating  lobe  would  form  for  an  equally  spaced  array  with 
spacing  equal  to  the  average  spacing  Dn  of  an  unequally  spaced  array;  that  is,  to  an  angle  £  ^ 
from  the  main  beam  given  by 

i  +  sin  J  =  jj-yx  •  (36> 

Beyond  this  region,  the  sidelobes  tend  to  be  more  or  less  random  with  an  average  relative  side- 
lobe  level  given  by  the  reciprocal  of  the  number  of  elements.  Thus,  control  over  a  sizeable 
fraction  of  visible  space  requires  that  the  average  spacing  between  the  elements  be  less  than  a 
wavelength.  For  such  spacing,  coupling  effects  are  certainly  not  negligible. 

A  practical  model  for  the  excitation  circuit  of  a  typical  array  of  independently  fed  elements 
was  indicated  in  Fig.  22.  Since  it  is  both  desirable  and  common  to  attempt  to  control  only  the 
open-circuit  voltages  of  the  generators  of  such  circuits  (rather  than  the  currents  which  depend 
on  the  mutual  coupling),  we  will  assume  that  these  voltages  have  been  chosen  to  be  equal  in  mag¬ 
nitude  and  progressively  phased  for  beam  pointing;  i.e.,  en  =  eQ  exp[-jkzn  sin{  ]  for  a  desired 
pointing  angle  {  ,  for  elements  located  at  points  zfi  along  a  line.  In  the  absence  of  coupling,  the 
array  factor  would  then  simply  be  proportional  to 

A(t-)  =  £-■  £  exp[jkzn(sin{  -sin4o))  (37) 

*  a  n 

We  have  seen  in  Fig.  41  that  the  agreement  between  the  array  factor  as  given  by  Eq.  (37)  and 
the  array  factor  as  computed  for  practical  arrays  including  coupling  is  quite  poor,  particularly 
in  the  sidelobe  region  near  the  main  beam.  The  reason  is  not  hard  to  find. 

As  pointed  out  in  Sec.  III-B-2,  the  contribution  of  an  element  to  the  far  field  at  any  angle  is 
proportional  to  the  pattern  of  that  element  at  that  angle  when  the  element  is  excited  in  the  pres¬ 
ence  of  the  remainder  of  the  elements,  the  latter  passively  terminated  in  Z^.  The  shape  of  such 
patterns  in  equally  spaced  arrays  depends  markedly  on  the  spacing  between  the  element  and  its 
neighbor.  Figure  43  shows  the  general  behavior  of  the  patterns  of  parallel  dipoles  0  25  x  above 
ground  in  a  linear  array  as  a  function  of  element  spacing  (the  behavior  of  other  elements  in  a 
linear  array  shows  a  similar  trend).  An  approximate  model  for  the  dependence  of  the  element 
gain  function  on  the  spacing  between  elements  of  a  linear  array ^  and  the  observation  angle  f  is 

l»t|lZ-<,[DA(f^|.lntl)l  MU'  <“> 

where  [x]  indicates  the  largest  integer  in  x,  and  |feU)|Z  is  the  power  pattern  of  an  isolated  ele¬ 
ment.  The  denominator  appears  complicated,  but  merely  represents  the  number  of  grating  lobes 


f  For  planer  array*,  the  variation  it  with  the  area  allotted  rite  element*  rather  than  with  the  length.  The  following 
argument*  are  applicable  to  such  array*  with  this  change  in  dependence. 
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Fig.  43.  Gain  functions  for  elements  in  linear  arrays  with  various  inter-element  spacings. 


ANGLE  (<*«> 

Fig.  44.  Comparison  of  parallel  dipole  array  patterns  obtained  using  both  exact 
and  approximate  methods  for  including  coupling  in  an  unequally  spaced  array. 
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in  visible  space  for  an  array  with  spacing  D  and  pointing  angle  4-  Figure  43  compares  this 
model  with  the  accurately  computed  element  patterns,  and  the  agreement  is  generally  good. 

This  agreement  suggests  that  the  patterns  of  unequally  spaced  arrays  including  coupling  may 
be  well  approximated,  if  the  change  in  spacing  is  a  slowly  varying  function  of  position,  by  a  sim¬ 
ilar  function  for  the  contribution  of  each  element.  Under  this  assumption,  the  array  factor  is 
given  by 


A(4)  *  2 


°n/X 


1 1/2 


1  7  [DnA(TT7"s'inlT)l)  exp[jkzn(sin{  -sin«0)] 


th 


where  is  the  spacing  allotted  the  n  element 


D  = 
n 


zn+l  zn-l 


(39) 


In  Fig.  44,  we  show  an  array  pattern  computed  with  the  use  of  the  approximation  of  Eq.  (39) 
as  compared  with  the  computed  parallel  dipole  array  pattern  of  Fig.  41.  The  comparison  with 
the  curve  of  Fig.  41  that  includes  mutual  coupling  effects  indicates  that  the  approximation  is 
quite  good.  A  further  simplification  is  suggested  by  the  observation  that  the  coupling  princi¬ 
pally  affects  the  sidelobes  near  the  main  beam  and  has  little  effect  on  the  far-out  lobes.  This 
suggests  that  substituting  the  pointing  angle  £q  in  place  of  4  in  the  denominator  of  Eq.  (39)  would 
be  an  acceptable  approximation;  i.e.. 


A<4)  *  Z 

n 


D  /A  1/2 

1  +  [DnAd  +'  I  sin  40f))  exp[jkzn(sin4  -sin40>] 


(40) 


Array  factors  computed  using  Eq.  (40)  compare  reasonably  well  with  those  computed  using  Eq.  (39). 


C.  Coupling  Effects  in  Arrays  with  Nonisolating  Feeds 

When  a  finite  array  is  fed  by  a  reactive  network,  several  effects  can  occur  in  addition  to  the 
gain  and  impedance  variations  discussed  above.  In  particular,  there  are  problems  associated 
with  multiple  reflections,  phase  shifter  properties,  resonance  phenomena,  and  bandwidth.  Some 
of  these  problems  will  be  discussed  in  this  section  for  both  reactive  corporate  feeds  and  reactive 
series  feeds.  The  corporate  feed  that  we  consider  is  constructed  with  three-port  (reactive)  junc¬ 
tions,  as  shown  schematically  in  Fig.  45  for  a  1:8  power  divider. 

J/ 

The  scattering  matrix  for  a  1:2  (=N)  reactive  corporate  feed  using  3-db  junctions  can  be 

obtained  by  following  a  wave  through  the  structure  when  only  one  port  is  excited  with  all  others 
properly  terminated.  Part  of  such  an  incident  wave  will  be  reflected  back  to  the  driven  port  by 
the  junctions,  while  the  remainder  will  continue  through  the  structure  to  the  other  ports.  The 
amplitude  and  phase  of  a  wave  out  of  a  particular  port  will  then  be  the  corresponding  scattering 
coefficient  when  the  input  wave  has  unit  amplitude  and  zero  phase.  For  example,  consider  a 
unit  wave  incident  on  an  output  port  of  the  1:8  divider:  at  the  first  junction,  half  the  incident 
wave  will  be  reflected  back  to  the  driven  port  and  will  arrive  with  (-2k i^-t)  phase  shift,  and 
similarly  for  the  other  reflections.  The  wave  returned  to  the  driven  port  for  this  case  can  be 
shown  to  be 

-j  exp(-j2kl1]  {1+7  exp(-j2ki2)  +  exp(-j2M*2  +  *3)|} 
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Fig.  45.  Model  for  a  1:8  reactive  corporate  feed. 

For  the  general  case  of  a  corporate  feed  with  K  levels  of  power  division  (1:2  power  divider), 
the  self-scattering  term  at  each  output  port  SQ  is 


S 


o 


-f  exp  l-j2kl  j] 


(41) 


The  scattering  coefficients  for  waves  coupled  from  the  driven  output  port  to  the  other  output  ports 
are 


S 


2n-i 


(42) 


where  2n  —  1  is  the  minimum  number  of  junctions  encountered  in  traversing  the  feed  from  the 
ith  output  port  to  the  jth  output  port.  (Note  that  n  ranges  from  one  to  K,  giving  only  K  distinct 
cross -coupling  terms.  Also,  for  n  =  K,  the  summation  term  in  brackets  is  zero.)  This  leaves 
just  the  forward  coupling  from  the  input  port  to  each  of  the  output  ports: 


S 


f 


=  2-K/2 


exp 


(43) 


It  is  readily  evident  from  the  above  results  that  the  scattering  matrix  for  this  type  of  feed 

is  strongly  dependent  on  the  electrical  lengths  of  the  junctions.  For  example,  the  self-scattering 

term  can  vary  in  amplitude  from  a  minimum  value  of  i/2  (for  f  ,  =  \/ 4  and  t  =  \/ 2.  n  >  2)  up 

to  a  maximum  value  of  1  —  l/Z1^  (for  i  =  X/2,  n  >  2).  The  other  scattering  coefficients  have 

K  ** 

similar  variations  ranging  from  1/2  up  to 

-  -J,  .  ,44, 

Also  observe  .hat 


lS2n-il  = 
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S1  ~  So  =  ^pHZWj] 


which  means  that  minimum  self-reflection  occurs  when  the  first-order  cross -coupling  is  max¬ 
imum  and  vice  versa.  These  limiting  values  for  the  scattering  coefficients  indicate  that  the 
scattering  matrix  can  vary  substantially  as  the  frequency  is  varied  over  a  moderate  range  when 
practical  line  lengths  are  used.  This  indicates  that  there  will  be  a  bandwidth  limitation  to  any 
practical  corporate  feed  system. 

With  the  feed  scattering  matrix  thus  defined,  a  relatively  straightforward  computation  will 
solve  for  the  overall  performance  of  the  antenna  plus  feed  network  system  for  a  specified  set  of 
parameters.  These  computations,  although  conceptually  quite  simple,  are  very  lengthy  even 
with  a  high-speed  digital  computer  because  inversions  of  two  large  matrices  are  involved  at 
each  scan  angle.  For  exact  calculations,  computer  storage  limits  us  to  relatively  small  arrays 
(«  64  elements);  however,  the  approximate  techniques  discussed  above  can  be  successfully  used 
to  solve  for  the  performance  of  a  few  elements  in  an  array  which  is  very  large  in  one  dimension. 

For  example,  we  have  examined  the  behavior  of  the  center  row  of  32  elements  in  a  32  x  65 
dipole  array  driven  by  65  identical,  uncoupled/  reactive  corporate  feeds.  When  reciprocal 
phase  shifters  are  used,  reflection  sidelobes  can  result  as  the  beam  is  steered  from  broadside, 
as  previously  discussed.  Figure  46  shows  an  example  of  these  multiple  reflection  effects  for  a 
square-grid  array  with  inter-element  spacings  of  0.6x  and  mounted  Q.25X  above  a  ground  plane. 
Note  that  the  reflection  lobes  in  this  case  are  very  high  because  of  relatively  large  antenna  mis¬ 
matches;  for  smaller  angles  of  scan,  the  mismatches  and,  hence,  the  reflection  lobes  are  de¬ 
creased.  If  nonreciprocal  phase  shifters  are  used  in  this  same  array,  the  reflection  lobes  are 
essentially  eliminated  as  indicated  in  Fig.  47. 

We  have  also  analyzed  an  8  x  8  square-grid  dipole  array  with  each  row  of  the  array  fed  by 
a  1:8  reactive  corporate  feed.  For  this  study,  the  array  had  an  inter-element  spacing  of  0.6X 
and  no  ground  plane  was  used.  An  array  of  this  size  exhibits  substantial  edge  effect  when  driven 
with  independent  sources;  however,  when  properly  designed  reactive  corporate  feeds  (t  f  =  l  2  = 

I  j  =  2.25X  in  this  case)  are  used,  the  currents  in  any  row  are  constrained  to  have  nearly  identical 
amplitudes  by  the  feed  network  as  demonstrated  by  the  top  curve  in  Fig.  48. 

This  is  a  direct  consequence  of  a  property  of  lossless  three-port  junctions  with  output  arm 
lengths  an  odd  multiple  of  a  quarter-wavelength  long:  the  amplitudes  of  the  currents  in  the  two 
finite,  nonzero  loads  terminating  the  output  lines  are  in  a  constant  ratio  (unity  for  3-db  junctions) 
independent  of  the  values  of  the  loads.  Unfortunately,  this  result  only  applies  for  a  discrete  set 
of  line  lengths  and/or  a  discrete  set  of  frequencies.  For  example,  a  5-percent  decrease  in  the 
operating  frequency  of  the  feed  described  above  results  in  the  array  illumination  shown  in  the 
bottom  curve  of  the  figure. 

When  this  array  is  scanned  by  nonreciprocal  phase  shifters*  to  point  a  beam  at  ©  =  45', 

<f>  =  45”,  the  current  distribution  across  half  the  array  at  f  =  0.95  f  is  as  shown  in  the  up¬ 
per  part  of  Fig.  49.  The  current  distribution  at  broadside  is  also  shown  in  the  figure,  from 
which  it  is  readily  seen  that  scanning  the  array  to  point  a  beam  off  broadside  can  substan¬ 
tially  alter  the  array  illumination.  On  the  other  hand,  if  the  array  is  operated  at  the  center 
frequency,  this  statement  does  not  apply,  as  demonstrated  by  the  dashed  curve  of  Fig.  50. 

t  That  it,  the  65  reactive  feed*  can  be  considered  a*  isolated  from  one  another,  as  would  be  the  case  if  each 
reactive  feed  were  fed  by  its  own  source  or  by  a  well-matched  hybrid  corporate  feed. 

t  The  phase  shifters  were  assumed  to  be  digital  with  a  smallest  phase  increment  of  45*. 
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Fig.  46.  Array  pattern  of  a  32-element  linear  array  fed  by  a  reactive  corporate  feed  network 
with  reciprocal  phase  shifters.  *• 


Fig.  47.  Array  pattern  of  a  32-elen.ent  linear  array  fed  by  a  reactive  corporate  feed  network 
with  nonreciprocal  phase  shifters. 
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Fig.  48.  Frequency  dependence  of  The  element 
currents  in  an  8x8  array  fed  by  eight  1:8  re¬ 
active  corporate  feed  networks. 


Fig.  49.  Scan  angle  dependence  of  the  element 
currents  in  an  8x8  array  fed  by  eight  1:8  re¬ 
active  corporate  feed  networks. 


Fig.  50.  Effect  of  feed  network  type  and  phase  shifter  type  on  the  element 
currents  in  on  8  X  8  array. 
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Fig.  51.  Cylindrical  beamforming  Ism  geometry. 
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Fig.  52.  Pattern  nomenclature  for  c  cylindrical 
beomforming  lent. 
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Fig.  53.  Change  in  geometry  with  feed  location*  for  a  cylindrical 
beamforming  lent. 


This  figure  also  shows  the  current  distributions  which  result  when  analog  reciprocal  phase  shift¬ 
ers  and  an  independent  (hybrid)  source  feed  are  used.  Note  that  the  use  of  reciprocal  phase 
shifters  results  in  a  very  poor  array  illumination  tape"  and,  as  previously  discussed,  can  cause 
very  high  reflection  sidelobes  when  the  beam  is  scan  led  off  broadside.  Also  observe  that  reac¬ 
tive  feed  networks  with  nonreciprocal  phase  shifters  can  give  an  m  ray  illumination  superior  to 
that  obtained  with  independent  sources  over  a  small  band  of  frequencies. 

IV.  EFFECTS  OF  MUTUAL  COUPLING  ON  RADIATION  PATTERNS 
OF  MULTIPLE-BEAM  OPTICAL- TYPE  ANTENNAS 

There  are  several  techniques  for  the  realization  of  simultaneous  multiple  beams  from  an¬ 
tennas,  with  perhaps  the  oldest  being  the  use  of  multiple  feeds  in  conjunction  with  an  optical- 
type  antenna,  such  as  a  reflector  or  a  lens.  For  example.  Fig.  51  illustrates  a  two-dimensional 
version  of  such  an  antenna  -  a  Luneburg  lens  fed  by  multiple  sources. 

The  multiple-feed  problem  has  a  precise  counterpart  in  the  theory  of  large  arrays,  and 
knowledge  about  the  effect  of  coupling  on  the  performance  of  large  linear  and  planar  arrays  can 
be  applied  directly  and  easily  to  this  problem.  To  illustrate,  let  us  analyze  the  two-dimensional 
(cylindrical)  Luneburg  lens  system  of  Fig.  51  (i.e.,  the  figure  is  assumed  to  extend  to  z  =  *•»). 

Consider  the  arrangement  of  Fig.  52.  The  lens  is  fed  by  a  single  feed  which  produces  an 
electric  field  per  unit  current  into  the  feed  of  e(v),  where  v  is  the  angle  measured  from  the  phase 
center  of  the  feed  (we  assume  it  has  one).  The  lens  performs  a  linear  transformation  on  e(y)  to 
produce  a  secondary  illumination  E(y).^  For  example,  if  the  cylindrical  lens  is  designed  so  that 
the  feed  phase  centers  are  on  the  lens  surface,  geometrical  optics  will  show  that 

E(y)  =  - ,  y  =  p  sin  v  .  (45) 

(pQ  cos  v)1/£  ° 

22 

The  far  field  at  a  range  R  is  then  determined  to  within  a  phase  factor  by 


E(y)  exp  (jky  sin*)  dy 


(46) 


where  k  =  2ir/\. 

The  actual  form  of  the  transformation  is  of  little  consequence  to  the  argument,  however. 

It  suffices  for  now  to  use  as  a  starting  point  the  statement  that  energizing  a  single,  isolated  feed 
at  a  location  corresponding  to  the  n  =  0  position  of  Fig.  51  produces  an  equiphase  illumination 
along  the  aperture  plane  of  E  (y)  per  unit  feed  current. 

If  the  single  feed  is  moved  to  a  new  position  -  say,  the  n  of  Fig.  51  -  the  aperture  plane 
moves  around  the  lens  through  an  angle  nD/pQ  radians,  as  illustrated  by  Fig.  53.  The  phase 
front  is  tilted  by  nD/pQ  radians  relative  to  the  y-axis,  and  the  phase  center  of  the  aperture  plane 
has  moved  from  x  =  Rq,  y  =  0,  to  x  =  pQ  cosnD/pQ,  y  =  pQ  sinnD/pQ.  Obviously,  if  nD/pQ  «1, 
the  displacement  of  the  phase  center  in  the  x-direction  becomes  negligible,  while  the  y  displace¬ 
ment  approaches  nD.  However,  if  pQ»  as  assumed,  most  of  the  pattern  structure  of  interest 
will  lie  near  {  =  0  of  Fig.  52,  and  the  y-offset  of  the  phase  center  will  have  little  effect  on  close- 
in  sidelobe  levels. 


t  The  dimensions  on  e(v)  and  E(y)  differ,  so  that  squaring  each  and  dividing  by  the  intrinsic  impedance  of  free 
space  produces  watts  per  unit  sol  id  angle  and  watts  per  unit  surface  area,  reqsectively. 
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Fig.  54.  Assumed  feed  circuit  configurations  for  a  cylindrical  beamforming  lens. 


As  a  result  of  the  foregoing  approximations,  the  change  in  the  aperture  illumination  on  a 
line  parallel  to  the  y-axis  that  results  from  moving  the  feed  appears  as  simply  a  change  in  phase 
front  slope 

En(y)  =  Eo(y)  exp  |jkn  Hz]  (47) 

when  viewed  at  small  angles  from  broadside. 

Returning  to  the  multiple-feed  configuration  of  Fig.  51,  if  we  excite  feed  0  with,  for  example, 
a  voltage  generator  of  open -circuit  voltage  V  as  in  P  ig.  54(a),  and  terminate  the  other  feeds  as 
in  Fig.  54(b),  then  the  current  in  the  n  n  feed  is  just 


I  =  V  Y 
n  o  no 


(48) 


where  Y  is  the  mutual  admittance  between  feeds  n  and  0  (a  function  of  the  feed  types,  feed 
no  r 

placements,  and  and  Z t  of  Fig.  54).  The  total  illumination  along  a  line  parallel  to  the  y-axis 
Et(y)  is  then 

N 

E,<y>  =  2  '„E„w 

n=-N 

which,  by  Eqs.  (47)  and  (48),  gives ^ 


N 


Et(y)  =  Eo(y)  VQ  £ 
n=-N 


(49) 


Thus  for  a  large  (pQ  »  optical  system,  the  far  field  in  the  region  near  the  main  beam  will  be 
the  Fourier  transform  of  Eq.  (49),  where  EQ(y)  is  the  aperture  illumination  per  unit  current  into 
the  feed  that  would  result  from  a  single  feed  placed  on  the  x-axis. 

Equation  (49)  can  be  put  explicitly  in  the  form  of  the  primary  feed  illumination.  The  effect 
of  the  optical  system  on  a  single  point  source  feed  is  to  transform  the  feed  primary  pattern  —  say, 
eQ(i>),  where  v  is  the  angle  from  the  x-axis  as  indicated  in  Fig.  52  —  into  EQ(y)  by  some  linear 
transformation,  which  we  will  denote  as  an  operation  I, 

E(y)  =  L  e(u)  .  (50) 

For  example,  for  a  surface-focused  Luneburg  lens,  the  transformation  is  given  by  Eq.  (45).  The 
exact  form  of  the  transformation  is  not  of  immediate  consequence  here,  only  the  fact  that  it  is 


tThe  purist  may  object  to  the  implicit  assumption  that  the  pattern  which  results  from  a  feed  being  driven  at  its 
terminals  and  the  same  feed  being  parasitically  excited  is  necessarily  the  same.  This  assumption,  which  was 
used  to  simplify  notation,  is  not  really  necessary  to  the  argument. 
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Fig.  55.  Spherical  beamforming  lens  geometry. 

linear,  which  allows  us  to  transform  both  sides  of  Eq.  (49)  to  give. 

N 

et(y)  =  eQ(1')  VQ  Tj  Yno  exp(jknD  sinr]  (51) 

n=-N 

where  we  have  used  y  =  p  sinr  in  the  exponent  as  appropriate  to  the  Luneburg  lens,  since  the 
phase  slope  in  Eqs.  (47)  and  (49)  was  derived  for  this  configuration. 

The  physical  interpretation  of  et(i>)  is  that  eQ(v)  of  Eq.  (51)  is  the  pattern  of  one  element  of 
a  linear  array  taken  in  the  presence  of  the  remainder  of  the  elements  when  they  are  passively 
terminated,  i.e.,  the  "element  pattern"  of  Sec.III-B,  if  there  are  enough  feeds  on  either  side  of 
the  element  in  question  to  assume  an  infinite  number.  The  results  are  directly  applicable  to 
three-dimensional  optical  feeds  as  in  Fig.  55,  with  the  planar  array  gain  function  used  for  feeds 
arrayed  as  in  the  figure. 

Thus  the  conclusion  of  the  analysis  is  that  a  multiple  feed  in  an  optical  system 
will  radiate  a  pattern  whose  shape  close  to  the  main  beam  is  given  by  trans¬ 
forming  the  pattern  of  one  feed  in  the  presence  of  the  other  feeds  when  they 
are  terminated  in  the  impedance  normally  used. 

This  conclusion  is  based  on  the  following  approximations: 

(a)  pQ  »  \,  and  therefore  geometrical  optics  is  valid. 

(b)  The  lens  is  large  enough  so  that  the  feed  assembly  (or  the  extent  of  the 
feed  over  which  elements  couple  appreciably)  is  small  enough  for  the 
curvature  to  be  negligible. 

(c)  The  array  is  large  enough  so  that  the  main  beam  and  sidelobe  region  of 
interest  is  of  small  angular  extent,  since  for  angles  very  far  from  the 
line  through  the  excited  feed  (the  x-axis  of  Fig.  51  if  feed  0  is  excited), 
the  difference  in  phase  center  locations  of  the  aperture  planes  for  the 
various  feeds  is  not  negligible  in  the  direction  of  interest. 

In  substance,  these  approximations  amount  to  a  restriction  to  very  large  (narrow-beam)  systems. 
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